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Abstract 


There are many attributes of electric propulsion which make it a desirable technology for 
use in modern warship designs. However, current motor technology makes electric propulsion 
noncompetitive from a cost standpoint. The single largest impediment to making electric 
propulsion more affordable is the size of the propulsion motors. Two technologies aimed at 
reducing the size of ship propulsion motors are investigated, rotor cooling via radial rotating 
thermosyphons and the multi-disk axial gap geometry. 


Size reductions in electric motors are ultimately limited by the ability to remove heat from 
the windings. Two-phase thermosyphons are considered for cooling the rotor windings of an 
axial gap motor because they can transfer large amounts of heat with a relatively small 
temperature difference. Predictions for the heat transfer coefficients of two-phase radially 
mounted thermosyphons are developed and experimentally evaluated. 


The multi-disk axial gap geometry significantly reduces both the weight and volume of an 
electric propulsion motor over conventional radial gap designs. This is accomplished by 
consolidating several machines together on a single shaft that share a common magnetic circuit. 
This novel geometry also allows significant reduction in the motor's diameter, giving the ship 
designer more flexibility in locating the propulsion motors within the ship. A computerized 
propulsion motor design tool is developed for conducting preliminary design studies of this type 
of motor. This tool can be used in conjunction with a ship synthesis model to conduct feasibility 
level ship studies in order to evaluate the total ship impact of this innovative motor technology. 


The two technologies described above have a synergistic effect when combined into a 
single design. A destroyer-sized motor of 35,000 HP @ 164 RPM exhibits a 15% weight 
reduction and more than a 30% volume reduction over the latest designs currently being 
developed by the U.S. Navy. 
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Title: Professor of Electrical Engineering 
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Chapter 1: Introduction 


1.1 Background and Motivation 


Due to the current political and budgetary climate, the United States Navy is under 
great pressure to reduce both the acquisition and operational cost of its ships. One 
technology which shows great promise for reducing the size and manpower requirements 
of ships, and consequently their cost, is integrated electric propulsion. Some of the 
advantages of this type of propulsion system include increased arrangement flexibility, 
increased fuel efficiency, increased automation and reduced manning, all of which 
contribute to reducing the cost of a ship without sacrificing its mission capability. 

The U.S. Navy has been conducting research in various technologies that will help 
make modern electric drive ships feasible in the near future. As part of that research, this 
thesis looks at novel cooling techniques as applied to axial gap motors that are very well 


suited to ship propulsion applications. 


1.1.1 Electric Drive Ships 


Electrically driven ships are not new. They have existed in one form or another for 
most of this century. However, since the 1940's electric drive has not been competitive 
with mechanical drive systems for warship applications. The variable speed requirement 
for ship propulsion systems has historically limited electric drive systems to the power 
levels achievable with brush-commutated direct current drive motors. Today this limit 1s 
only about 10,000 horsepower. A typical 8000 ton destroyer requires 30,000 - 40,000 


horsepower per propeller shaft. Modern computerized control systems and high power 
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semiconductor switching devices have made variable speed electric drive systems possible 
using wound-field or permanent magnet synchronous motor technology. 

The decision process of what goes where on a modern warship is largely a 
trade-off between many conflicting design requirements. Propulsion and electrical 
systems, while essential for every ship, are not part of the mission payload and thus their 
location within the ship should be subordinate to the placement of the payload items. In 
the case of a warship, those payload items are the combat systems. 

Mechanical drive systems must maintain a stringent alignment between the 
propulsion prime movers and the ship's propellers. This constraint severely constricts ship 
arrangement options. Additionally, it makes construction more difficult and expensive due 
to the long propulsion shafts which must be aligned during construction. 

Steam driven ships require boilers, condensers, steam turbines and other very large 
and heavy propulsion equipment. To ensure the stability of the ship, these heavy items are 
placed low in the ship's hull. This arrangement works well since this location 1s required 
anyway to maintain the mechanical drive system's alignment. In the past thirty years, the 
U.S. Navy and most other navies have switched to gas turbine mechanical drive on 
combatant ships. Gas turbine engines are very small and light when compared to a steam 
propulsion system of similar power rating, but they have other drawbacks. Gas turbines 
require large amounts of air, are non-reversible and suffer from poor fuel efficiency when 
operated at off-design conditions. 

In modern mechanical drive ships, the gas turbines are located low in the hull to 


maintain alignment with the propeller shafts (see Figure 1-1a). This arrangement results in 
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Figure 1-la: Conventional ship arrangement. 





-Propulsor ship arrangement. 


Figure l-1b: Podded 
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long intake and exhaust ducts which take up much valuable space aboard the ship. The 
non-reversibility of the engines is usually compensated for via reversible pitch propellers. 
While reversible pitch propellers have worked very well in the fleet, they are costly and are 
inherently less efficient than a properly designed fixed-pitch propeller. 

All of these drawbacks could be eliminated with an electric drive system. By 
eliminating the mechanical connection between the gas turbines and the propellers, the 
engines could be placed higher in the ship where they would require less volume for intake 
and exhaust ducting. This results in a smaller ship design. The variable speed and 
reversing requirement can be handled via the control of the semiconductor switching 
devices which feed power to the propulsion motors. This eliminates the need for 
reversible pitch propellers. Combining the propulsion and ship's service electric power 
loads on a common prime mover reduces the total number of installed gas turbines on the 
ship. If they are sized properly for the ship's operating profile a significant fuel savings can 
be expected. Electric drive also provides the ability to "cross-connect" the propulsion 
motors; that is, to split the output from one operating gas turbine generator between both 
propulsion motors in a twin screw ship. This flexibility existed with steam propulsion and 
was lost upon switching to mechanical gas turbine drive. 

With present day technology, electric propulsion is still not cost competitive with 
modern mechanical drive systems. This is primarily due to the large size and weight of a 
direct-drive electric motor. In order to achieve all of the benefits cited above, the 
propulsion motors would also have to be moved outside the ship's hull as depicted in 


figure 1-1b. This concept is known as "podded-propulsors" and further restricts the size 
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of the propulsion motor beyond that required for an "inside-the-hull" electric drive 


afrangement. 


1.1.2 Propulsion Motors 


Electric motors for ship propulsion are subject to a rather unusual set of 
constraints. They can be best characterized as low-speed high-torque machines. Recall 


that power (P) and torque (T) are related by: 


P=T-Q (ied) 


where 2 1s the rotational speed in radians / second. A typical propulsion motor for an 


8000 ton destroyer might be rated at 30 MW @ 164 rpm, a torque of 1.75 x 10° N-m. 


Torque in a conventional radial gap electric machine is given by: 


eer Lt (1.2) 


where T is the electromagnetic shear stress, R and L are the rotor radius and length, 


respectively. In an axial gap machine, the torque is given by: 
_2_(p3_p3 i 
T=2n(R3-R})-1 Gi) 


where R, and Rare the outer and inner radii of the rotor. Equations (1.2) and (1.3) 
assume one air-gap in the motor. In either geometry, the radius should be as large as 
possible to maximize the generated torque. Normally, the limit to the radius is determined 
from structural considerations. However, the ship's geometry constrains the radius to less 
than structural limits would normally allow. For inside-the-hull use, the outer radius of the 


motor would be limited to about 2 meters. For a podded-propulsor ship, the outer radius 


ie? 








could be no more than | meter with a length limit of about 2.0 meters. Current motor 
technology can meet the in-hull limit; however, the size and weight make the system 
noncompetitive with a mechanical drive system from a cost standpoint. The size 
limitations of the podded-propulsor configurations are beyond current motor technology. 

There has been a great deal of effort aimed at reducing the size of propulsion 
motors. Numerous technologies have been suggested for use including superconducting 
homopolar DC and water-cooled permanent magnet AC machines [1-1],[1-2]and [1-3]. 
Common to all of the proposed technologies is the need to use sophisticated cooling 
techniques. 

As can be seen from (1.2) and (1.3), to reduce the size of any electric machine 


without reducing the torque output, the electromagnetic shear stress must be increased. 


The electromagnetic shear stress (t) is given by: 


t=B,-Kg (1.4) 


where B, is the air-gap magnetic flux density, and K, is the stator surface current density. 
The magnetic flux density is limited by saturation effects in the iron of the machine, so the 
only way to further increase the electromagnetic shear stress is to increase the current 
density. Increasing the current density also increases the resistive heating losses in the 
windings. Therefore, the current density is limited by the ability to remove the generated 
heat from the windings. For a ship propulsion motor to be successful, it must have high 
current densities and consequently requires some sort of advanced cooling methods. 

The cooling methods which are typically used in large machines are not suitable for 


shipboard use. The most common methods used today in large turbine generators are 
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hydrogen cooling and direct water cooling. Due to the explosive nature of hydrogen gas, 
it is unacceptable for shipboard use. Helium could be considered as a substitute; however, 
its density/specific heat product is only 2/3 that of hydrogen gas. This means 50% more 
volumetric flow rate is required to remove the same amount of heat with the same 
temperature difference. 

Water cooling could be used and has been considered in previous propulsion 
motor designs. In direct water cooling the conductors of the winding are hollow and 
cooling water passes through these hollow conductors in direct contact with the copper 
conductors. The same cooling water must also come into contact with pumps and heat 
exchangers which are at ground potential. The main drawback to direct water cooling is 
the high water purity requirements which must be maintained for the water to span the 
required voltage potential without shorting out the windings. While technically feasible, 
direct water cooling poses the potential to become a major maintenance problem at a time 
when the Navy is looking to significantly reduce its maintenance expenses. 

Another cooling method exists which has become the central focus of this 
research. Two-phase rotating thermosyphons are able to move heat away from the rotor 
winding to the central part of the rotor nearly isothermally. Heat fluxes on the order of 
10° Watts/meter’ are possible with these devices. The heat generated in the windings may 
then be removed from the motor by passing cooling air or water through the center of the 
motor. Two-phase capillary heat pipes could be used in a similar fashion to remove heat 


from the stator windings to the periphery of the machine. Although this complicates the 


design and construction of the motor, the thermosyphons and heat pipes are virtually 
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maintenance free once installed. As a result, this type of motor would be much simpler to 


operate and maintain than one which requires direct water cooling. 
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Figure 1-2: Multi-disk axial motor geometry. 


Another technology which shows promise for ship propulsion applications is the 
multi-disk axial gap motor. Figure 1-2 shows a schematic representation of this geometry. 
With the radius constrained the only way to increase the size of a conventional machine 1s 
to lengthen it. With the outer radius similarly constrained, the size of an axial machine 
may be increased by reducing the inner radius or by combining multiple machines on the 
same shaft. It turns out that the optimum ratio of r,/r, for maximum torque generation 1s 
in the range of 0.5-0.6 depending upon the number of pole pairs [1-4]. The multi-disk 
geometry is simply the consolidation of several axial gap machines together on the same 


shaft that share a common magnetic circuit. Although this does increase the length of the 
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machine, this increased length is much less than that required for a conventional radial gap 
machine to achieve a similar power rating. 

The combination of the multi-disk axial gap geometry with the highly efficient 
cooling capability of radial rotating thermosyphons and heat pipes will allow significant 
size reductions in the ship propulsion motor. These size reductions are necessary to make 


electric drive cost competitive with mechanical systems. 


1.1.3 Thermosyphons and Heat Pipes 


Thermosyphons and heat pipes have been in existence for many years. Both are 
highly efficient heat transfer devices which operate by recirculating a fluid within a closed 
hollow tube. Figure 1-3 shows a schematic representation of the thermosyphon. Heat is 
transferred via latent heat during the evaporation and condensation of the working fluid. 
While a thermosyphon relies on gravity to return the liquid portion of the working fluid to 
the evaporator, a heat pipe contains a wick structure that uses capillary action to return 
the fluid to the evaporator. This allows a heat pipe to work in zero gravity or inverted 
(with the condenser section below the evaporator section). Both thermosyphons and heat 
pipes can achieve an effective thermal conductivity which is several hundred times greater 
than that of copper [1-5]. 

For analysis purposes, the thermosyphon can be separated into its constituent 
components: condenser, evaporator and adiabatic sections. Vapor enters the condenser 
section where it is condensed on the walls via the film condensation process. The liquid 


film flows through the adiabatic section into the evaporator on the walls of the heat pipe. 
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When the liquid reaches the evaporator section, the wall temperature rises above the 
saturation temperature and the film begins to evaporate. At some point in the evaporator, 
the film will reach the pool of liquid in the bottom of the evaporator, replenishing it. The 
liquid pool undergoes a pool boiling process, creating vapor to repeat the cycle. Although 
treated as separate phenomena in most heat transfer texts, in the thermosyphon each 


process 1s related to the others through the conservation laws. 


Condenser 
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ee | Section 
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Figure 1-3: Typical thermosyphon. 


The radial rotating thermosyphon (RRT) is related to the "rotating heat pipe" 
which was first introduced by Gray [1-6] in 1969. Both devices rely on the centrifugal 
force of rotation to return the condensate to the evaporator. In the rotating heat pipe the 


axis of rotation is parallel to or coincident with the pipe axis. In the radial rotating 
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thermosyphon the axis of rotation is perpendicular to the axis of the pipe. This simple 
geometry change significantly alters the analysis of fluid flow and heat transfer in the 
device. 

Capillary and rotating heat pipes have been developed for many applications and 
their heat transfer characteristics are fairly well understood [1-5]. Radial rotating 
thermosyphons have been largely ignored to date. One of the contributions of this thesis is 


to provide some insight into the heat transfer characteristics of these devices. 


12 Previous Research 


There are several fields of research which apply to the present study. These 
include ship propulsion, axial gap motors, cooling of electric motors, heat pipes & 
thermosyphons, two-phase flow, heat transfer and gas turbine technology just to name a 
few. The following section summarizes some of the more important and interesting works 
in the various fields of applicability to this thesis. 

There have been numerous papers published concerning the design of axial gap 
motors. The vast majority of these deal with small permanent-magnet A.C. and D.C. 
machines. Takeda [1-7] discusses a variable-reluctance machine and Varga [1-8] has 
published several nearly identical papers about induction machines. All reviewed 
publications on axial gap machines concentrate on the basic electromagnetic circuit and 
design performance of the motors. Little attention is paid to the structural or thermal 
design of the motors in these papers. None of the papers concerned with axial gap motors 
consider anything other than air-gap cooling. Di-Napoli [1-4] develops relationships for 


torque, iron & copper weight, iron & copper losses, power and efficiency vs. radius ratio. 
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Chan [1-9] and Desequilles [1-10] both discuss the possibility of using a multi-disk 
geometry with permanent magnet machines. Chan also calculates the optimum radius 
ratio for maximum power and minimum rotor inertia. Several of these authors hint at or 
state that axial gap machines are more power dense than conventional machines, but none 
prove this assertion. The paper by Leung & Chan [1-11] is the only paper reviewed that 
considers space harmonics of the magnetic field in the air-gap. This paper discusses pole 
shaping on salient pole synchronous axial machines. 

There have been numerous publications concerning thermal modeling of 
conventional radial gap motors: [1-12], [1-13] and [1-14] for example. However, there 
seem to be no published studies of the thermal behavior of axial gap motors. The 
underlying principles are the same, but the geometry is significantly different which makes 
any models developed for radial gap machines invalid. 

A rotating disk that 1s convectively cooled on its sides 1s the basic geometry to be 
considered for modeling the air-gap heat transfer in the axial gap machine. There has been 
an abundance of experimental and theoretical work done on this problem, both with and 
without a nearby stator. Most of the research was concerned with the cooling of gas 
turbine rotors. The authority on this particular heat transfer problem appears to be J.M. 
Owen who has published several papers on the subject. The most useful of his works is 
[1-15] which provides a complete set of design formulae for determining heat transfer 
coefficients for both the stator and rotor as a function of rotational speed, air-gap length 


and air mass flow rate. 
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There have been a few papers published regarding electric machines with heat-pipe 
or thermosyphon cooling: [1-16], [1-17], [1-18], [1-19], and [1-20]. All of these are 
exclusively concemed with radial gap machines. Consequently the heat pipes considered 
are axially oriented on the rotor of the machine and are only of limited use in the present 
study. 

The field of literature which has been useful in the present study of radially 
oriented thermosyphons is gas turbine blade cooling. In the late 1940's and 1950's there 
was quite a bit of research done on liquid and two-phase cooling of gas turbine blades. 
Developments in this area apparently died out as high temperature alloys were developed 
for turbine blades. There was a resurgence of interest in this area in the late 1970's largely 
due to the energy crisis. Of the early work, the most interesting paper is by Cohen & 
Bayley [1-21] that considers using radially mounted two-phase thermosyphons to cool gas 
turbine blades. They built a rotating test rig to show that they could transfer heat from the 
tip to the center of the rotor, but all of their quantitative experiments were carried out on a 
fixed gravitational thermosyphon. The later work was concentrated on the more 
fundamental problems of flow and heat transfer in rotating fluid films. Of particular 
interest here are the theses by Mudawar [1-22] and El-Masri [1-23]. Both of these works 
look at heat transfer from thin rotating films. 

There is an extensive body of literature on thermosyphons and heat pipes. 
Although there are numerous references in this area, most of the pertinent research can be 
summarized with just a few papers. The paper by Al-Farah [1-24] is the most complicated 


numerical analysis of a tilted thermosyphon. It models the thermosyphon using a 
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three-dimensional finite difference solution. Harley [1-25] provides a fairly useful 
transient two-dimensional numerical model for thermosyphons. Reed & Tien [1-26] 
develop a very useful one-dimensional analytical model for the thermosyphon which 
includes a prediction of the flooding failure mode. Shiraishi [1-27] does one of the better 
experimental analyses and develops empirical equations for the heat transfer coefficients 
based upon Nusselt theory. The paper by Nguyen-Chi [1-28] contains a good description 
of the various operating limits for thermosyphons. This paper also develops the flooding 
limit theoretically and corroborates it experimentally. 

Most of the other papers are quite repetitive, with a modified Nusselt analysis 
compared to their experimental results. All of the previously published work in the 
thermosyphon field has been done in a gravitational system. Apparently, there has been no 
published research concerning thermosyphons in a centrifugal gravity field. 

Another related field which has seen much press Is that of the "rotating heat pipe" 
which is actually an axially mounted rotating thermosyphon. P.J. Marto has published 
numerous papers on this subject, [1-29] provides a very good overview of the research 
through 1984. Because of the geometry difference between the axial and radial rotating 
thermosyphons, the research in the rotating heat pipe field cannot be applied to the radial 


rotating thermosyphon. 


1.3 Aim of Present Research 


The research underlying this thesis hopes to accomplish three main objectives. The 
first objective of this thesis is to develop a computerized design tool for conducting a 


preliminary design of multi-disk axial gap motors. This design tool will include 
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electromagnetic, thermal and structural aspects of the motor design. In order to properly 
develop the thermal design of the motor, it is necessary to consider both the air-gap 
cooling and cooling of the windings via advanced cooling methods such as thermosyphons 
and heat pipes. Herein lie the other two objectives of this research. 

An overall thermal model of the axial gap motor based on classical heat transfer 
models will be developed for use in the computerized design tool. Such a model currently 
does not exist for axial gap motors. Development of such a model is one contribution of 
this research. Part of the overall thermal model will be the thermal modeling of the radial 
rotating thermosyphon for rotor cooling. Since there has been no prior research into the 
thermal behavior of the radial rotating thermosyphon, the thermal model which 1s 
developed will be evaluated experimentally. This development of the theory for radial 


rotating thermosyphons will be another contribution of this research. 
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Chapter 2: Axial Propulsion Motors 


Axial gap motors represent a significant departure from conventional radial gap 
electric motors. Figure 2-1 shows a comparison of the geometry. The primary advantage 
of this geometry in ship propulsion applications is the ability to increase the active volume 
of the machine by adding disks which is much more volume efficient than adding length to 
a conventional radial gap machine. This chapter develops a model of the multi-disk axial 


gap machine for preliminary design purposes. 


Stator 


4 _ Back Iron 


_— — Rotor =... 
Winding r 


_—>— Stator Winding 


Structure 


Axial Motor 





Figure 2-1: Radial Vs Axial Motor Geometry 
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2.1 Voltage and Current 


For the preliminary design analysis of the axial gap machine, only the space 
fundamental component of the air gap magnetic flux will be considered. The air gap is 
assumed to be small such that the air gap flux exists in the axial direction only. Under 


these assumptions, the air gap magnetic flux is given by: 
B(a) = Bmaxsin p(@ — ot) (2x) 


The flux per pole 1s obtained by integrating equation (2.1) over the active area of the 


machine: 


n/p °O 
= [°° |” B(@)r- dr do. (2.2) 
the result is: 


— (12 - 1?) BmaxC0S (pot) (2.3) 


The terminal voltage can be obtained by recalling that: 


v= aie where A= ®-N, - Ss, or in the steady state: Vins = ee 
/2p 


dt 
where N_ is the number of series turns per stator phase and S, is the number of stators 


connected in series. The terminal voltage 1s then: 


ON,Sskokp( 12 - 7?) Bee 


Vims = (2.4) 
/2p 


where the breadth and pitch factors, k, and k,, have been introduced to account for 
distribution and chording of the windings. The breadth and pitch factors are calculated in 
the same manner as for a conventional radial gap machine. See Fitzgerald, et. al [2-1] for 


breadth and pitch factor calculation methods. 
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The phase current 1s obtained next in terms of the stator current density. From 


geometric considerations, the current is given by: 


aa Netots : A slot ; 


I, 
6-N, 


Ja-Sp (2.5) 


where S, is the number of stators connected in parallel. The stator space factor (A,) is 


now defined as the fraction of the stator circumference which contains conductor slots: 


Nstots "Ws 
= 2.6 
: OTT. ( ) 


where w, 1s the width of a stator slot. Note that this parameter has been defined in terms 
of the inner radius. This avoids the problem of having the teeth go to a zero or negative 
thickness which could occur if 4, is defined at the mean or outer radius. The area of one 


slot is given by: 
Astlot = hs "Ws (2.7) 


Equations (2.6) and (2.7) are then substituted into equation (2.5) to obtain the phase 


current: 


a artiAsh, ace (2.8) 


where the product (h, J,) is simply the electric loading, K,. The terminal variables of the 
machine are now defined by equations (2.4) and (2.8) in terms of the electric loading, 
magnetic loading and geometric parameters. The power rating of a three-phase machine 1s 


then: 


P =3V,-I,-cos(y) (2.9) 
where wis the power factor angle. 
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2.2 Reactance and Resistance 


Next, the reactance and resistance values for the windings must be determined. 
Only two reactances are considered for preliminary motor design, the synchronous 
reactance and the leakage reactance due to slot leakage. End-turn leakage, zig-zag 
leakage and differential leakage reactances are left to a more detailed design analysis. The 
slot leakage reactance can be calculated by integrating Ampere's Law around the contour 


shown in Figure 2-2: 

f H-di= { 7 en (2.10) 
which becomes: 

Hew), - w-h, rl) 
The total flux crossing the slot is simply the product of the flux density and the area: 

®) = WoH- hg (ro — 13) (2.12) 


The leakage reactance is defined by: 


oN,D 


Xal = i. 


(uley 


Equations (2.8), (2.11) and (2.12) are substituted into equation (2.13) to obtain the 


required expression for the leakage reactance: 


7 6N2kpk2 10h, (ro salir 


= 2.14 
F 2nrjAsSp ea 
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Figure 2-2: Ampere's Law Contour for Slot Leakage Calculation 


The synchronous reactance can be obtained by considering the flux linked by the 
stator winding due to its own excitation. The MMF produced by exciting a single stator 


phase 1s: 


r= 7 sin(p8) (2.15) 


* 2p 
The corresponding flux density in the air gap 1s: 


— 2ptoNsla . 


y= re Doe (p9) (2.16) 


The flux per pole is found by integrating from the inner radius to the outer radius and over 


the angle of one pole pitch: 
n/p (To 
b= |)" [°B.-rdr do (2.17) 


Carrying out the integration yields: 
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2woNgla(r? ~ r?) 


®D = 2.18 
ae (2.18) 
which gives the phase reactance as: 
2wooN2k?2k?(r2 - 3) 
X .49 = ————_— (2.19). 
MEP 

The synchronous reactance for a three-phase machine is determined by: 

Xa = 2Xaa0 + Xa (2.20) 


2 


where X,,, and X,, are given by equations (2.14) and (2.19). The synchronous reactance is 


per-unitized by dividing equation (2.20) by the base impedance: 
Be Nib 
Ip 


where V, and I, are the base voltage and current, respectively. 


Ly (22) 


The required field excitation for the motor 1s determined from the phasor diagram 
at rated conditions, see Figure 2-3. Assuming the base voltage and current are the rated 


values, the per unit field excitation is: 


ear = [1 +x2 +2xqsin(w) (2) 


The field current required to develop rated voltage at no load 1s: 


7B maxPS (2.23) 


Ie = 
= 2uokrNe 


where k, is the field winding factor and N, is the number of series turns in the field winding 


referred to the stator. The rated field current is then: 


Tp = Tin + Car (2.24) 
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Figure 2-3: Phasor Diagram at Rated Conditions 


Finally, the field current density is given by: 


Nel¢ 
i ei >) 


where A, is the rotor space factor corresponding to A, on the stator and h, is the rotor slot 
height. 

The turn length must be determined in order to calculate the winding resistances. 
Figure 2-4 shows the typical winding geometry for an axial gap machine. As a Starting 
point, the end turns are assumed to be triangular with approximately 45 degree bend 
angles. Although this is not exact, it is a good approximation for a motor with a large 


number of poles as are being considered. This gives the turn length as: 


Li = 2(ro—1i) + ¥2 (to +1aL (2.26) 


ue 
p 


Equation (2.26) holds for both the stator and rotor windings where "pitch" is the 


where « is the angle between corresponding sides of a coil and is given by: a = =- pitch. 
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fractional pitch (i.e.. 5/6 pitch) which is unity for the field winding. The resistance for one 


turn of a coil is: 


a 
Figure 2-4: End Turn Geometry (Armature Winding Shown) 


Li 


vo 
Gai Awire 


(227) 


where Ocu is the conductivity of copper [mho/m]. A,_. is the cross-sectional area of one 


nue 


turn of the winding and is given by: 


2nrjh,A,SSF 
wire = —— 2.28 
‘i 2N,:m ( ) 


where m is the number of phases and SSF ts the fraction of the slot cross-section which 
contains conductors. The stator space factor (SSF) and slot height are replaced with their 


respective rotor values for the field winding. The winding resistance is then found from: 
R=R,N;,S, (2.29) 
Combining equations (2.27) - (2.29) yields: 


mL, -N?2-S, 


= es 2.30). 
neha SSE cey) 
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Equation (2.30) is valid for either stator or rotor windings when the appropriate values are 


used. 
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Figure 2-5: Section View of Motor 


2.3 Structural Design 


For preliminary design purposes, the structural design of the motor will consider 
only the major components of the machine. Specifically, the shaft, rotor disks, stator disks 
and case are sized based upon classical engineering mechanics. While an analytical 
analysis such as this is rather simplistic in nature, it serves the purpose of generating 
reasonable dimensions for comparison purposes in the preliminary design stage. Figures 


2-5 and 2-6 depict the geometry under consideration. 
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Figure 2-6: Detail of Stator and Rotor Construction 


The shaft is modeled as a hollow steel tube subjected to torsional loading. Popov 


[2-2] provides a simple design formula based upon purely elastic behavior: 


I a “ 
a. ~ Tmax (2.31) 








where d, is the outer diameter of the shaft, T is the load torque, Tmax 1S the maximum 


allowable shear stress of the material and I, is the polar moment of inertia which is given 


by: 
Ip = (ds -d8) (2.32) 


for a hollow cylinder. Torque is equal to P/Q, so if we assume that d = 0.8d, then the 


equations (2.31) and (2.32) can be combined and solved for d, as: 
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Pp 1/3 
= — 4 


For design purposes, tmax for steel is assumed to be 100 MPa. 


The stator and rotor disks must be strong enough to transmit the electromagnetic 
torque to the motor case and shaft, respectively. The required area for transmitting the 


torque from the disk to the drum or case 1s found from: 


T 
——— 


ar Tmax 


(2.34) 


where r 1s the inner radius for rotor disks and the outer radius for stator disks and T,, 1s the 
torque per disk. The rotor is constructed from solid steel, so t_,. for the rotor will be that 

of steel. The structural portion of the stator disk is constructed from a composite material 
whose maximum shear stress is about 14 MPa. T,, is determined from the rated power, 


speed and number of disks as: 





Tee (2.35). 


The required disk thickness 1s then determined from: 


eee (2.36) 
2 
For the stator, r will be: 
ee stator pitch (2.37) 
2 
to account for the end turns. Similarly, for the rotor r 1s given by: 
pe rotor pitch (2.38) 


2 
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In the region of the rotor disks, the shaft outer diameter is increased to that of the 
inner diameter of the rotor disks, creating effectively a "drum" region to support the rotor 
disks. This is done both to reduce rotor weight and to allow more space for cooling 
devices in the inner region of the machine. Equations (2.31) and (2.32) can be solved for 


the drum inner diameter in terms of the known outer diameter. The result is: 


, 16dor T 1/4 
Gir = ds = oT (2.39) 


where T_,, 1S the maximum shear stress for steel and equation (2.38) determines the radius 
corresponding to d,,. 

The case that supports the stator disks is sized in a similar manner and ts also 
constructed of steel. Unfortunately, d, is now known and d, 1s the unknown in 
equations (2.31) and (2.32). The inner diameter is twice the radius given in equation 


(2.37). The result is an implicit equation in terms of the outer diameter: 


d5s— Fare doe = di (2.40), 


An implicit equation is not the preferred result, but it is easily handled in the design 
spreadsheet. This is discussed further in Section 5.1. 

This simple structural analysis doesn't begin to address all the concerns of 
designing a motor, but it does allow rationally based sizing of the major structural 
components of the motor. Obviously, a more complete structural analysis which considers 


deflection of the various components, stress concentration due to geometry changes and a 
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dynamic analysis to determine shaft critical speeds should be performed before proceeding 


further with any motor design. 


2.4 Weight and Volume 


An estimate of the weight and volume of a particular design must be obtained in 
order to make comparisons of various machine designs. For preliminary design purposes, 
only the electromagnetic and structural components of the motor are considered in these 
estimates. 


The stator winding volume is calculated from: 


and the rotor winding volume is similar with the appropriate rotor values substituted for 
the stator parameters. The windings are imbedded into slots in magnetic iron on both the 
rotor and stator disks. The volume of the stator magnetic material can be approximated 
by: 

vai = (12-17 )(1—As)Ns - (2h, +s) (2.42) 


where A, is the stator space factor defined by equation (2.6), h, is the stator slot height 
and t. is the disk structural thickness determined from equation (2.36). Since the density 
of the structural composite portion of the stator 1s small, its volume 1s not calculated 


explicitly. The stator case volume is given by: 


Vease = E42, ea) -L (2.43) 
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where d., 1s the case outer diameter calculated in equation (2.40), d.. is twice r., from 


equation (2.37) and L is the stack length of equation (2.49). 


The rotor volume 1s calculated as the sum of the disks and the drum. This volume 


Vi =n (13,-13) te (12-17) 2h, |-N, — Vreu + (42, — a?) (14) 


where r,, and r,, are the structural radii determined from equations (2.37) and (2.38), 
respectively. d,, is simply twice r.. and equation (2.39) determines d... 

There must be back iron at both ends of the machine to provide a return path for 
the magnetic flux. The thickness of the back iron is determined from saturation limits in 


the iron and the desired flux density in the air gap. In the air gap, the flux per pole 1s: 


(7 a a Bag 


2 ae 4 
© ; (2.45) 


where B,, is the flux density in the air gap. In the back iron, the flux per pole is: 
® = (ro — Ti) Batt (2.46) 


which must be the same flux as that given in equation (2.45). Equating these fluxes and 


solving for the back iron thickness results in: 


Bag i at 
meee (OT) 2.47). 
tb Bex D ( ) 


The back-iron thickness is inversely proportional to the number of pole pairs. The 


back-iron volume is then: 


Vbi = an(r? om 2) tb (2.48). 
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Another important parameter is the stack length of the axial machine. The stack 


length (L) 1s: 
L=N,(t, +2h,)+Ns(ts + 2h;) + 2th + Nog (2.49). 


In all designs currently being considered, N,, N, and N, will be equal. 
The weight of the machine is simply the summation of the various volumes 


multiplied by the appropriate density of the component materials. 


2.5 Losses and Efficiency 


There are several loss mechanisms in electric motors. The primary ones are the 
winding loss, core loss, windage and friction loss. These three mechanisms are considered 
in this analysis. Winding loss is simply the I°R loss associated with the current passing 


through the armature and field windings. This is calculated as: 
Py = 312Ra+1;Re (2.50) 


for a three-phase machine. Winding loss is the single largest loss component. This is 
fortunate as this is the only loss mechanism which can be calculated in a straightforward 
manner. 

Core loss arises from eddy currents induced in the magnetic materials of the 
machine by the magnetic fields created in the windings. The specific core loss is a function 
of flux density in the core, electrical frequency, lamination thickness and resistivity of the 
core material. Walker [2-3] provides empirical graphical relationships for core loss as a 
function of lamination thickness, electrical frequency, steel grade and flux density in the 
core. For the typical values of 1.3 T flux density at 60 Hz in 0.5 mm laminations of grade 


170 steel, the specific core loss is approximately 33 KW/m’ from his graphs. The total 
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core loss for a particular machine is the product of this specific core loss and the volume 
of the core material. To account for hysteresis effects, the total core loss is multiplied by 
an empirical factor of 1.55. 

Windage and friction losses arise from the motion of the rotor. Air in the air gap 
and the bearings which support the rotor both cause a drag on the rotor as it rotates. 
However, due to the very slow speed of the propulsion motors being considered, this loss 
is very small in comparison to the other loss mechanisms described above and can be 
ignored. 

The efficiency (n) of the machine is then given as the ratio of the power output to 


the power input. This can be wnitten as: 


T] = ] _ 2 Prosses es 1) 
Pin 
where ZP,..... is the summation of the winding and core losses. P,, is the power input to 


the motor terminals. 


2.6 Heat Transfer 


Heat transfer to the surroundings is the method by which all losses are ultimately 
dissipated. There are several paths through which heat may be transferred out of the 
machine, all of which act in parallel. Some of these paths occur naturally such as 
conduction out of the shaft or natural convection off the outer casing. The most 
important paths are the ones designed into the machine specifically for the purpose of 
cooling the heat generating parts. This cooling is provided via the air gap and any liquid 


or gas cooling passages designed into the machine to cool the windings. 
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2.6.1 Air Gap Heat Transfer 


Virtually all electric machines push air (or some fluid) through the air gap to 
provide cooling to the machine. In this respect, the axial motor is no different. However, 
the geometry of the air gap 1s significantly different from a conventional radial gap 


machine. This section examines heat transfer in the axial geometry air gap. 
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Figure 2-7: Axial Gap geometry 


The geometry of the axial gap motor is depicted in Figure 2-7. This heat transfer 
problem has been studied in gas turbine rotor cooling applications. Owen [2-4] provides a 
review of the theoretical and experimental work done in this area. The heat transfer in the 
rotor-stator system depends primarily on the rotation rate, gap ratio and the amount of 


throughflow. The rotation rate is measured via the Rotational Reynolds Number: 


or 
Vv 





Reg = 
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where @ is the rotational speed, r, is the outer disk radius and v is the kinematic viscosity. 
The gap ratio (G) is defined as the ratio of the axial spacing between the rotor and stator 


and the outer radius of the rotor: 
_&§ 
St 
where g is the gap length. The amount of throughflow, Q is given by: 
Q = 2nrieVs 


and is typically non-dimensionalized as the flow coefficient, C,: 


Q 


rv 


Cw 





The throughflow can also be non-dimensionalized as the radial Reynolds number, Re:: 


for comparison with the rotational flow. D, is the hydraulic diameter which is defined as 


four times the cross-sectional flow area divided by the wetted perimeter: 








4A. 4-2nrg 
D =o = = 2 
7 Pp 4nr 5 


and is twice the gap length for the geometry under consideration. 
For heat transfer between the rotor and the air, the heat transfer coefficient is 


defined non-dimensionally by the Nusselt number: 





0.8 
Nu, = ia 2 0.145(S) (2.52). 
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where h, is the dimensional heat transfer coefficient for the rotor with SI units of 
[Watts/(meter -°C)] and k is the thermal conductivity of the fluid in the air gap. The total 


heat transfer from the rotor is: 

Q=h,-A-(T;- Tin) (2.53) 
where A is the surface area of the rotor, T, is the rotor temperature and T._ is the inlet air 
temperature. Similarly, the heat transfer coefficient for the stator is defined by: 


Niele 


Nus = - 





= 0.76 - C257. G-032 (2.54). 


The total heat transfer between the stator and air gap is calculated using (2.53) when the 
appropriate stator values are used. Equations (2.52) and (2.54) were developed by Owen 


[2-4]. 


2.6.2 Direct Water Cooling 


Direct water cooling is one of the most effective winding cooling methods used 
today. As the name implies, the cooling water comes into direct contact with the windings 
of the machine. This is accomplished by using hollow conductors and pumping the 
cooling water through them. Because the heat generated in the windings does not have to 
be conducted through electrical insulation before reaching the cooling water, this method 
iS very attractive for high voltage machines which have thick winding insulation. The 
drawback to this method is that the cooling water must span the potential difference 
between ground and the winding voltage. This requires very high purity water, which 


from a shipboard maintenance standpoint is a serious deficiency. 
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2.6.3 Indirect Water Cooling 


Indirect water cooling is similar to direct water cooling except the cooling water 
stays outside the electrical insulation. The advantage of this method is that the high purity 
water requirements of direct water cooling are no longer required. The disadvantage is 
that heat must be transferred through the electrical insulation before reaching the cooling 
water. For a given maximum winding temperature, this reduces the allowable temperature 
rise of the cooling water. The motor must be designed with larger cooling passages and 


lower current densities which increases the machine size. 


2.6.4 Hydrogen Cooling 


While water cooling ts very effective for stator cooling, the hydrostatic pressure 
caused by the centrifugal force of rotation makes it impractical for rotor cooling in most 
applications. Hydrogen, due to its low density and very high specific heat, 1s very well 
suited to rotor cooling. Hydrogen is passed through hollow conductors in a direct cooling 
mode or it can be used in lieu of air in the air gap and indirect cooling passages built into 
the machine. This cooling method is common in large turbine generators. The drawback 
of this method is the flammability and explosive nature of hydrogen gas. This potential 


safety hazard makes hydrogen cooling unacceptable for shipboard use. 


2.6.5 Air Cooling 


The advantage of air cooling is that it is simple and reliable. Lower heat transfer 


coefficients are obtained than with other methods. Current densities are necessarily 
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limited to much lower values, resulting in a much larger and heavier machine. Air-cooled 


motors do not provide the high power densities needed in ship propulsion applications. 


2.6.6 Two-Phase Cooling Methods 


All of the winding cooling schemes previously described have relied upon 
transferring heat from the windings to a cooling medium via sensible heat, or increasing 
the temperature of the coolant. Two-phase cooling relies on the heat of vaporization of 
the coolant to remove heat from the windings. There are two advantages of this cooling 
method. First, the large heat transfer coefficients obtained via boiling of the coolant allow 
higher heat fluxes into the coolant. Second, the need for a temperature difference in the 
coolant is virtually eliminated, allowing the designer to use more of the total allowable 
temperature difference for conduction through the electrical insulation to the coolant. 

While two-phase cooling has been previously used to cool conventional machines, 
it has yet to be considered for cooling axial gap machines. The designs developed for this 
thesis use two-phase rotating thermosyphons to cool the rotor windings. Development of 


the necessary design equations 1s deferred to Chapter 3. 
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Chapter 3: Radial Rotating Thermosyphons 


There are essentially three heat transfer processes which occur within the radial 
rotating thermosyphon (RRT). These are film condensation, film evaporation and pool 
boiling. This chapter develops a description for these processes which is applicable to the 


RRT and also investigates the limits to thermosyphon performance as applied to the RRT. 


3.1 Film Condensation 


The classical heat transfer problem of film condensation of a vapor on a vertical 
wall is known as a Nusselt analysis and is treated in any introductory heat transfer text 
such as Mills [3.1]. The primary assumptions built into the Nusselt analysis are that the 
inertia of the film is neglected and the temperature profile is linear. The effects of vapor 
superheat, liquid subcooling and vapor velocity are also neglected. While seemingly quite 
restrictive, the standard Nusselt type analysis gives good results in many applications. In 
the RRT, gravitational force is replaced by the centrifugal force of rotation which may be 
One or two orders of magnitude larger and varies with radius. If the cross-sectional area 1s 
very small, the vapor velocity may become significant at higher heat transfer rates. With 
these differences in mind, this section develops a Nusselt type analysis that 1s valid for the 
variable centrifugal force found in the RRT. The momentum and energy equations are 
then evaluated to determine the range of validity of this solution. 

Analysis of film condensation begins by looking at an elemental volume of fluid 
within the film (see Figure 3-1). A force balance per unit width on this element in the 


radial or x-direction yields: 
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Figure 3-1: Elemental Film Volume 
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After rearranging and taking limits, this becomes: 


Sa + (01~ Pats) = 0 (3.1) 


This differential equation for the velocity profile 1s subject to the boundary conditions: 


by 


u(y = 0) = 0 (no slip at wall) 


and os =0 @ y=65 (neglect interfacial shear stress). 


Equation (3.1) can be directly integrated to yield the film velocity profile: 


u(x,y) = cor inan Oo -y?) (3.2) 


The mass flow rate of the film per unit width (I) is obtained by integrating equation (3.2) 


across the film thickness: 
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Energy Balance 





Figure 3-2: Energy & Mass Balance on Film Strip 


5 Xx =r ey 6? 
P= J} piu(x,y)-dy = SOP theim-s (3.3) 

The film Reynolds number (Re) is defined as: 
Ree (3.4) 


= 
Equations (3.3) and (3.4) can be re-arranged to give the film thickness (5) in terms of the 


film Reynolds number: 


344;V, Re aS 
SS oro) 
4(p1 — oa | a 


An energy balance on the film element of Figure 3-1 results in the differential 


equation: 


OT 


ay? = 0 (3.6) 
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which 1s subject to the boundary conditions: 


T(y = 9) = Taat 
T(y = 0) =Tw. 


Integrating equation (3.6) results in the linear temperature profile: 


T- ls = (Ts a a) 
which indicates that heat is simply conducted through the liquid film. The heat flow 


through the film is given by: 


d= ki Elo = h(Tan - Tw) (3.7) 


where the local condensation heat transfer coefficient has been defined in terms of the heat 
flux and the temperature difference across the film. A mass balance on a thin strip of the 


liquid film as depicted in Figure 3-2 shows: 
5 WW 5 
J, pu(x,y)-dylx +h - Ax -[ piu(xy)- dylssax =0 


this can be solved for the condensation rate: 


; i edu 5 _ dv 
m = dx I, piu(x, y)- dy = dx (3.8). 


An energy balance on the same element of Figure 3-2 shows: 


5 WW 5 T 
bs p1u(x, yh) -dyl,x +m hy - Ax -| p(x, yh - dyl eax ki ly-0 = (3). 


If subcooling of the liquid film is neglected and a constant specific heat 1s assumed, then: 


OT dV 
ae oe) 


Now, equations (3.4), (3.7) and (3.9) can be combined: 


oD 








| 
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ee 


Rene a= heey | = (3.10). 


For the radial rotating thermosyphon, the Nusselt number is defined by: 


h J v/@ 


ky 





Nu = (3.11) 


where /v/@ is taken as the characteristic length scale. Equation (3.10) can now be 





written as: 
dx__ Pr dRe 
{vio ~ 4Ja Nu ols) 
where the Prandtl number and Jakob number have been defined as: Pr = a and 


= Co(T sat a ix) 





Ja 7 , respectively. An energy balance on the entire film from x = 0 to 
fg 
x = L, shows: 
1m Lihg, 
hL<(T sat a Tw) = ie ; hfe = 4 Kei, 


which can now be written 





and the average Nusselt number is defined as: Nu = 


h/vi/o 
ky 


in terms of the exit Reynolds number: 





— oe 3.13). 
Nu ie Aya Rei. ( 3) 


The difference between the gravitational and centrifugal forces can be accounted 
for using the method presented by Dhir & Lienhard [3.2] which uses an effective gravity. 


Using their method, the film thickness is given by: 
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1/4 
A4uiki(T sat = yx 


SS , (3.14) 
(oi Pv)pihp, Ser 
4/3 
- 5 = —XtgR) 
where: get = Fe R@ax (3.15). 
The local Nusselt number is given by: 
/ 1/4 
ge ae Pl Peg (3.16). 


=o tome 
k ° Aviki(T sat — Tw) 
For the case of the RRT, g(x) = ra? = X.@? + x@? where X, is the distance from the axis 
of rotation to the top of the condenser and x is the distance down from the top of the 
condenser. The radius of curvature (R(x)) is infinite for a flat surface. The resulting 


effective gravity is: 


2 
Ler = a OLE (3. l 7) 
]1-[X./(X, +x)] 


and for p; >> py, the local Nusselt number becomes: 


_[Xetx) (3 i 3.18 
) {Re (3.18). 


The average Nusselt number is obtained from equation (3.13). First, the exit Reynolds 
number must be found by substituting equation (3.18) into (3.12) and integrating along the 


length of the condenser. 


1/3 
i, (X_+x)'? -dx= a ( (vile ; ve fo“ (3Re) -dRe 


which integrates out to: 
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3/4 
4| (Xe+Le)"-X2 3ia 


3 ( vie ) Pr 


eT. = (3.19) 


Equation (3.19) is substituted into equation (3.13) to obtain the average Nusselt number: 





(3.20). 


/ = 
we( 2)" [Oterto”-xP]""_ A 
3Ja [ee ky 


Equation (3.20) can be used to estimate the heat transfer coefficient in the condenser of 
the radial rotating thermosyphon. A Nusselt number based on condenser length can also 


be defined as: 


Pr) 1/4 Oey eee 


a 4/3 


(3.21). 


3Ja 





An important effect that has not been included in this result is the Coriolis force. 
This is treated later. Subcooling, which has also been neglected can be accounted for by 


replacing the latent heat with a corrected latent heat given by: 
hy, = hg + 0.68C p(T sat - Tw). 


The above analysis is limited to laminar flow which occurs when the film Reynolds number 
is less than 30. 

For the wavy laminar regime in gravitational systems, the empirical relationship of 
Chun & Seban [3.3] is commonly used to predict the heat transfer coefficient. Their 


relationship is given by: 
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Figure 3-3: Evaporating film 


/ heg/vy 
Nu = —— = 0.822Re°” (3.22). 
l 





For use in rotational systems, the length scale is changed to /vi/@ . The resulting local 


Nusselt number 1s: 


( 1/3 
Nu = 0.822 Ream | ee | (3.23) 


J vi/o 


Equation (3.23) is then substituted into equation (3.12) and integrated over the length of 
the condenser to obtain the condenser exit Reynolds number: 


G 0.82 
3Ja_ [Xe +L) -xe | 


Pr ( [vie ) 4/3 


Rei, (3.24). 


The average Nusselt number is then found from equation (3.13) to be: 
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Nu Z 3( Pr 0.18 vile Cs Lene -x#2 


a\4Ja L. Tm (3725): 
( [vile ) 
Similarly, a Nusselt number based upon the condenser length is defined as: 
Pr.) O18) (Xe tLe) - XV 
= =| — ———— .26). 
AJa Cae 


aa 4/3 





This result is valid over the wavy laminar regime of 30 < Re < 1000-1800. The radial 
rotating thermosyphon is not expected to operate with a turbulent film so that regime will 
not be considered here. For turbulent films, the reader is referred to the works of 


El-Masri [3.4], Mudawar [3.5], Dakin [3.6] and [3.7]. 


3.2 Film Evaporation 


Film evaporation is treated in a manner similar to film condensation. However, the 
film thickness is not zero at x = 0. Ifthe adiabatic section is assumed to have no heat or 
mass transfer, then the condenser exit film Reynolds number will equal the evaporator 
entrance Reynolds number. Because the liquid film is evaporating, the energy balance on 


the film element shown in Figure 3-3 results in an equation similar to equation (3.12): 





vile 4Ja Nu 


dx — Pr dRe G2Ty 


Equation (3.18) is still valid here. It is substituted into (3.27) and integrated from x = 0 to 


x = L, and Re = Re,, to Re = Re,, yielding: 
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_4|(3 - — 3Ja (Ke+L.)° -x# 
en = 3 (3rex. Dr a (3.28) 
4g vi/@ ) 
The average Nusselt number is then found from equation (3.13): 
3/4 
Nea (Fe) . " Sia Cea) ee 
Nu in a qReLe : (3 29). 


Pr (Juve ) 4/3 


Just as in the condenser, another Nusselt number based on the evaporator length is defined 


as. 


4/3 _ y4l3 
_ 3Ja (Xe +e) — Xe" (3.30) 


Pr ( Jvio ) 4/3 





Equation (3.30) is used to predict the heat transfer from the film above the liquid 
pool in the laminar regime. For the wavy laminar regime, equation (3.23) represents the 
local Nusselt number vice equation (3.18). When this is integrated along the evaporator 


length the exit Reynolds number becomes: 


0.82 
122 3Ja_ (X-+L.)"°- X35 
Rete = | Re; 5, ine (3,51). 
( v\/@ ) 
Similarly, the average Nusselt number 1s: 
0.82 
43 yr43 
No = 21m Pr} peize _ 3a (Xe tbe) — Xoo (3.32) 


L. 4Ja = Pr ( [vio ) 
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for the wavy laminar regime. The Nusselt number based on the evaporator length is: 


4/3 +7 4/3 
elt = 33a (XetLe)  — Xe" (3.33). 


Pr ( [vio ) 4/3 





3.3. Nusselt Analysis Restrictions 

As stated at the beginning of this chapter, there are many assumptions built into 
the Nusselt analysis of Section 3.2. The conservation equations for a falling film in a 
rotational gravity field are now examined more closely. This is to determine when the 
analysis of Section 3.2 may be expected to accurately predict the heat transfer capability of 
the radial rotating thermosyphon. 

Figure 3-4 shows a thin film in a rotational gravitational field. A mass balance on 


the element in that figure gives the following continuity equation: 


Ou , OV 


xt eno (3.34). 


A momentum balance on the same film element taking into account both the w’r and the 


Coriolis-hydrostatic forces yields the following momentum equation: 


yt 4 yOu = 9 Babi Pv) 0? - 2oy (3.35). 


Ox Oy x oe 
An energy balance on the same element gives the energy equation as: 


2 
ee ee ge Gail (3.36). 


Ox Oy  y? 
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Figure 3-4: Film in Rotating Gravity Field 


These equations are nondimensionalized by defining the following nondimensional 








variables: 
é u * V * xX € y T ed Jie * ® * r 
u"=—: vam x*=—: = ; © = ———_; oO* = ne oe 
Uc? Ve’ lee ne lv,/o Tsat — Iw Od max Ro 


where the x-direction characteristic velocity is chosen to be the average film velocity: 


ue =U= 
p16 
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The y-direction characteristic velocity (v.) is the one unknown left to be determined. The 


continuity equation (3.34) 1s used to define v.: 





{vila 


Uc Ou" | Ve Ov" 














Ue du" ee. (3.37). 
lenox Milo Cy" EZ 
The momentum and energy equations can now be nondimensionalized as: 
. Lae TY Bye) = Bye Re One | 
2 
Reply * at GO ) = mleo 0) (3.39) 
Oy* Pry oy”? 


where the following non-dimensional groups have been defined: 


_ _ 4T _ Inertia Force 
Film Reynolds No. Re =  eonG aE SEes 


Be eis py) _ Buoyant Force 


Archemedes No. = oe ionar oes 
- _ _v1__ Viscous Force 

Berns. SaaS QL? Coriolis Force 

Beer Gl -p, _ Oe tiv; _ Momentum diffusivity 





k, ©! Thermal diffusivity ~ 


The order of magnitude for the above parameters 1s: 
10'< Re <10° 
10" < Ar, < 10" 
10°< Pr, <10' 


107< Ek, <10*. 
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Within this range of parameters, the viscous term and the gravity term in equation (3.38) 
have an order of magnitude of about one. The convective terms have an order of 
magnitude of 107 or less and can be neglected without significant error. The Coriolis term 
in equation (3.38) can also be neglected as it has an order of magnitude of 107 - 10°. 

In the energy equation (3.39), the Reynolds-Prantl number product must be small 
compared to 10° for the Nusselt assumptions to be valid. This condition is met for the 
fluids and operating conditions under consideration. Recall that the operating conditions 
which generate higher Prantl numbers are at low heat fluxes where the Reynolds number is 
at the smaller end of its range. 

At higher Reynolds numbers the Nusselt assumptions are no longer valid and the 
previous analysis can not be used. However, the gravitational term becomes smaller with 
respect to the other terms as the Reynolds number is increased. If this term can be 
neglected from the momentum equation then equations (3.38) and (3.39) resemble the flat 
plate boundary layer problem which is treated in many introductory heat transfer texts. 
Previous researchers [3.8] have applied the integral method for solving the classical 


boundary layer problem to gravitational thermosyphons. 


3.4 Pool Boiling 


Any portion of the liquid charge which is not in vapor form or forming part of the 
liquid film remains at the bottom of the evaporator forming the liquid pool. Heat transfer 
from the walls surrounding the pool is by natural convection at low heat fluxes and by 
pool boiling at high heat fluxes. Classical pool boiling can be described by the relation of 


Rohsenow [3.9]: 
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1/2 
where L, = | —s—_ , C,, and m are empirical constants that are a function of the 
Pi-Pv)gs 


liquid and surface characteristics. The relationship in Equation (3.40) cannot be used 
directly because the characteristic length is a function of gravity which varies with radius. 
Equation (3.36) can be modified by introducing a new characteristic length which is the 
same as that used in the condensation and evaporation Nusselt numbers above. The 


modified Nusselt number is: 


VO yy’ = Avo Jal ee et Pv 2202) Me 





Nu= 3.4] 
L. C3,Pr™ ( ) 
which is now a function of pool depth. The average Nusselt number is found from: 
=== J V1/o Ja’ & = ee M2 @ [lp 1/2 
i a (pee) ex (3.42) 
C3.Pr™ 2) Ly J ; 
which 1s: 
vi/@ Ja’ (p;-py\ !? (eps) Xe" 
2 Avila We (pi=pe)" | Op tba)" Xp? ce 
3 Cr Pr™ 21 Lp 


Equation (3.39) accounts for the effect of the variable gravity field in the liquid pool. A 


Nusselt number based on the evaporator length can be defined as: 





(3.44). 
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Figure 3-5: Body Force Components in Thermosyphons 


3.5 Coriolis Effects 


The Coriolis force occurs in any rotating system when the mass under 
consideration moves in the radial direction. A complete treatment of Coriolis acceleration 
can be found in [3.10]. In the case of the RRT the working fluid flows radially as the pipe 
is spun around, resulting in a Coriolis force as shown in Figure 3-5. The body forces of a 
tilted gravitational thermosyphon are also depicted in Figure 3-5. Lock [3.11], suggests a 
comparison of these forces as a method for accounting for the Coriolis force in the RRT. 

In the tilted gravitational thermosyphon the gravity force components on the liquid 


film are given by: 


2a =P1-g-cos(O) - gravity component parallel to thermosyphon axis 
2: =p1:g-sin(@) - gravity component perpendicular to thermosyphon axis. 
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where 0 is the angle between the thermosyphon axis and vertical. In the radial rotating 


thermosyphon the body force components are: 


ga(x) = (Pp1 — Pzg)@?(X +x)- centrifugal force 
21(x) = 2-p)-@- u(x) - Coriolis force. 


An effective tilt angle for the rotating thermosyphon can be defined as: 


2p i@u(x) 


8.) = ——__7_ 
re (P1—- Pg)@?(X +x) 


(3.45). 


where u(x) is the fluid velocity in the radial direction. Lock determines characteristic 
velocity scales from the equation of motion neglecting the convective terms. The resulting 
Coriolis tilt angle from his analysis is: 


1/3 
a: 13 D Ja Pi 
6. a (27) a Pr (p1 = —) (3.46) 


where small angles have been assumed. 

Lock further postulates that the Coriolis effects can be treated by looking at tilted 
gravitational thermosyphons. There have been several investigations of tilted gravitational 
thermosyphons [3.8], [3.12], [3.4] and [3.13], all of which show an increase in heat 
transfer ability with tilt angle up to a maximum that occurs in the range of 40-60 degrees 
from vertical. This increase is due to the thinning of the film layer by the gravity 
component normal to the thermosyphon axis. The film is thickened at the "bottom" of the 
circumference, but studies have shown this area to be very small in comparison to the 
thinned area; hence the resultant increase in overall heat transfer coefficient. For the 


conditions encountered in the RRT for axial motors, the Coriolis tilt will only be on the 


68 











order of one or two degrees using this method. This is consistent with the small size of 
the Coriolis term in equation (3.38). Such a small tilt angle will not have a significant 
impact on the heat transfer capability of the RRT. The asymmetry which it causes may 
help promote dryout failure on the leading edge of the evaporator. 
3.6 Performance Limits 

Dunn & Reay [3.14] provide a fairly complete accounting of performance limits 
for heat pipes. The performance limits of the radial rotating thermosyphon are similar to 
those of the heat pipe when the geometry differences and effects of rotation are accounted 
for. This section describes each type of limit and suggests a possible method of 


determining the limits for radial rotating thermosyphons. 


3.6.1 Viscous Limit 


The viscous limit of a thermosyphon occurs at low temperatures when the viscous 
drag dominates the vapor flow in the thermosyphon. This limit is directly proportional to 
the vapor pressure and density in the evaporator. A two-dimensional analysis of the vapor 


flow in the pipe yields the following result: 


2 
_ TyhgePvpy (3.47) 
1Opy Leg 


which agrees well with published data for heat pipes [3.11]. This limit is normally only 
encountered under conditions of very low vapor pressure and is not reached under the 


Operating conditions being considered. 
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3.6.2 Sonic Limit 


The vapor flow in a heat pipe or thermosyphon behaves much like the flow 
through a converging-diverging nozzle. The difference being that in the nozzle the mass 
flow is constant while the cross-sectional area varies; in the thermosyphon the area is 
constant while the mass flow rate varies. The evaporator exit is the point where the 
greatest vapor mass flow exists. This location corresponds to the throat of a nozzle. The 
sonic limit occurs when the vapor velocity at the exit of the evaporator reaches the speed 


of sound at that point. This choking limitation can be predicted by: 


q=0.474hg /pypy (3.48) 


which is cited in [3.11]. 


3.6.3 Entrainment Limit 


Entrainment is caused by the interfacial shear stress between the liquid film on the 
walls and the vapor core. The entrainment limit is often the limiting factor in low 
temperature thermosyphons such as are the subject of this research. Unfortunately, there 
seems to be no concise method for calculating this limit. Dunn & Reay [3.11] present the 


following correlation for thermosyphons: 


q=hefifofs (Pv [g(Pi-pwo) (3.49) 


where f, - f, are empirical constants which are functions of Bond Number, pressure 
coefficient and inclination angle, respectively. This relationship 1s used to predict 


entrainment failure in rotating thermosyphons. 
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3.6.4 Burnout Limit 


The burnout limit is encountered at high radial heat fluxes. Although this limit 
could be encountered in thermosyphons, it is primarily a problem in heat pipes when 
boiling within the wick structure upsets the capillary pumping action. In the RRT, this 
limit manifests itself in two distinct ways. One is the boiling crisis situation (or transition 
to film boiling) in the liquid pool at the end of the evaporator. The other is the critical 
heat flux limit in the liquid film above the pool. 

The geometry of the pool within the RRT is markedly different from that of classic 
pool boiling. Previous researchers [3.15] and [3.16] have shown that the shape of the 
pool boiling curve is preserved, but the critical heat flux is altered. Work with axial 
rotating heat pipes [3.17] has also shown that the centrifugal force can affect the critical 
heat flux. The critical heat flux can be calculated using the method of Lienhard and Dhir 


[3.5] for finite bodies. For the geometry of the RRT the critical heat flux is given by: 
max = 0.90 - Qmaxs = 0.90 - 7 py *hre[og(x)(P1 — oie (3.50) 


where q_.., is the flat plate pool boiling limit. By substituting in the centrifugal 
acceleration (g(x) =r?) for the gravitational acceleration (g), q,,,, has become a function 
of pool depth. 

Recalling that r= X, +x in the liquid pool, equation (3.50) can be integrated to 


obtain an average evaporator critical heat flux. The result is: 


5/ 
qmax = 0.90: Sythe [ow?(p1- pv)}! | (Xp + Lp) —X}* | (3.51) 


which accounts for the differences in both geometry and body force. Unlike most other 


limits on the RRT, the critical boiling heat flux is a radial flux value which acts over the 


WM 


entire heated surface area of the liquid pool rather than the cross-sectional area. As a 
result, other performance limits will be encountered well prior to the burnout limit of 
equation (3.51). 

The second type of burnout limitation takes place in the liquid film above the pool. 
Mudawar [3.5] studied evaporation and boiling of rotating liquid films in the field of gas 
turbine cooling. The critical heat flux in the liquid film occurs when the film boils away 
and can no longer re-wet the surface creating a dry wall condition. Based upon an 
analysis which balances vapor drag and body forces, Mudawar determined the critical heat 
flux in the liquid film to be: 


(PI —poeo i | (3.52) 


Vv 


which he confirmed experimentally. 
Another related phenomenon which is important to the radial rotating 
thermosyphon is the incipient boiling limit in the liquid film. This is the level of heat flux 


at which the liquid film begins to boil. This heat flux can be determined as: 


hee ky 


ee (T= Tex) 3.53 
SoT wav t) (358) 


i 


where v,, is the difference in liquid and vapor specific volumes. Rohsenow [3.18] bases 
this result on the premise that when the wall temperature exceeds a certain value, 
wall-bound bubbles will begin to grow. Since most commercially prepared surfaces have a 
wide range of cavity sizes, the incipient boiling condition should be independent of surface 


condition. 


a2 





Mudawar [3.15], presents equation (3.53) in a slightly different form: 


SS 3.54 
q ne (3.54) 
where the heat transfer coefficient, h is calculated from: 

h=0.042-A-Re®!7. Pr®3> (3755) 
an \ = Pr-0.48 Rel2!8(Pr-0.3)°!8-0.17] (3.56) 


The heat transfer coefficient in equation (3.55) is his theoretical result for turbulent films 
which is valid for 5 x 10°< Re < 80x 10°. For the conditions to be encountered in the 
present research, equation (3.54) can be used in conjunction with equations (3.29) and 


(3.32) to predict the incipient boiling condition. 


3.6.5 Dry-Out Limit 


A limit peculiar to thermosyphons 1s the dry-out limit. Dry-out occurs when the 
liquid film evaporates entirely prior to reaching the liquid pool at the bottom of the 
evaporator. It is caused by small fill ratios and long evaporator sections where the total 
heat transfer ability of the evaporator exceeds that of the condenser section. Dry-out can 
be predicted by comparing the instantaneous mass flow rates of the condenser and 
evaporator sections. When the mass flow rate of the evaporator exceeds that of the 


condenser, dry-out 1s imminent. 
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3.6.6 Flooding Limit 


Flooding or hold up 1s specific to thermosyphons. It is another result of the shear 
stress at the liquid/vapor interface. Flooding occurs when the vapor velocity increases to 
a point where the vapor drag on the liquid film overcomes the gravitational force, halting 
the downward flow of the film. Flooding differs from entrainment in that the film is not 
carried off in the vapor stream but 1s simply held stationary on the tube wall. This 
phenomenon 1s also common in two-phase vertical flow. Wallis [3.19] suggests the 


following relationship for predicting the heat flow at which hold-up will occur: 


De hte / 3Pv(P1 — Pv) 


Qinax = oy (6557) 
4 [1 +(pv/pi) "| 


When SI units are used, the empirical constant, C, = 1.105. Nguyen-Chi and Groll [3.20] 
achieved good results with this correlation on their tilted gravitational thermosyphons. 
Equation (3.57) is used to predict flooding failure in rotating thermosyphons with g 


replaced by the effective gravity due to the centrifugal force. 


3.7 External Heat Transfer 


Thus far, this chapter has dealt with the heat transfer processes within the 
thermosyphon. The thermosyphon must form part of a total cooling system to be of any 
practical use. This section examines how heat is transferred into the evaporator and out of 
the condenser of the thermosyphon. 

Heat is generated in the rotor both within the windings and in the magnetic iron. A 
two-dimensional finite element model of the rotor slot region was developed to determine 


the temperature distribution in the region of the rotor windings and thermosyphon 
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Figure 3-6: Rotor Slot Detail for Finite Element Solution 


evaporators. The geometry under consideration is depicted in Figure 3-6. This geometry 
simulates the rotor slot and surrounding iron at the inner radius of the rotor. The specific 
dimensions chosen are from the sample motor depicted in Figures 5-1 and 5-2. 

The temperature profile of Figure 3-6 indicates that the largest temperature drop 
occurs across the electrical insulation, all other temperature drops are small by 
comparison. This specific geometry results in a maximum winding temperature of 159° C 
when the wall temperature of the thermosyphon is assumed to be 100° C. 

At the opposite (condenser) end of the thermosyphon the heat 1s dissipated to a 
secondary cooling fluid. This secondary coolant may be air or water. Air provides fora 


simpler design, but more heat can be extracted with water than air. To determine if air can 


oi 
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be used to cool the thermosyphons the external heat transfer coefficient is calculated based 


on the classical flow through a tube bundle problem. 


Note: Conductor 
slots omitted for 


clarity. 


Rotor Disk 


View looking in radial direction of 


motor at average condenser radius. 





Figure 3-7: Condenser Section External Arrangement 


Figure 3-7 shows the geometry under consideration. The RRT's are assumed to 
have round cross-sections. The transverse spacing is taken as the average over the length 
of the condenser since it varies over the condenser length. The transverse and longitudinal 


pitch are calculated as: 


NM 
NM 


T pe 


Pr=5> 


L 
3 (3.58) 


The Nusselt number for the array of tubes is simply: 
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Nup 


Nup = ON (3.59) 


where N is the number of tubes (thermosyphons) Nup is the Nusselt number for a single 
tube and © is: 


CemeSe ons 
ee eee 3.60 
EES SID (CIT walle 


where: 


fat and . = 1 = ———- P, <1. 
eli ASP 





The single tube Nusselt number ts given by: 


0.62-Rey’Pri3 — Os sain 
D 


—— 
Nup = 0.3 +————2> —_. 
[1 +(0.4/Pr)?? |" 





— 0.62 - Rep? Pr" “ 
Nup = 0.3 + Eom | (_Rep _| 


+ -2x104 < Rep < 4x10° 
[1+(0.4/Pr)??] 282, 000/ | 





I 0.62 -Rel?pr!3 | (Rens), 


4/5 
Nios ‘4x10°<R 5x10°. 
"=D Te 4/Pr)?3 |" \282, 000/ | eu ee 


The Reynolds number is based on the tube diameter: 


where the velocity is the maximum between the tubes: 


pe: St 


= = mD ~ 
Seas 
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The above set of equations 1s used to calculate the maximum heat which can be transferred 
from the condenser to the air flowing through the motor without cooling fins. 

In most cases this heat transfer will be insufficient to completely cool the 
condenser. There are two possible solutions to this problem: 1) switch to water cooling 
of the condensers, and 2) add fins to the exterior of the condensers. If solution 1 1s 
chosen, then the above analysis still holds with water properties substituted for air in the 
equations. If solution 2 is chosen, then the problem is solved in a different manner. 


Figure 3-8 depicts the geometry of the condensers with cooling fins attached. 


View from end of motor: 
only 2 thermosyphons 
shown for clarity 


ee ee 


PTT > Thermosyphon 


REAR EE Condensers 
——— . Fins 


Enlarged View 


Tn ~ Centerline of 


thermosyphon 


Perimeter, P = 28+2(r2-r1) 
Wetted Area, A = S*(r2-r1) 


Fin Thickness, 2t 


eal Fin Spacing, S 


Figure 3-8: Finned Condenser Arrangement 
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Heat transfer from the finned surface is determined from the fin efficiency. For 
design purposes, the fins are assumed to be circular with an outer radius equal to half the 


distance between condenser tubes. The fin efficiency for such an arrangement 1s: 


nee iB iB 2) Gy 2) (6.61) 
: a) Ko(B-1r1)-11(B-r2) + 1o(B +1) - Ki (B- r2) | 


where I and K are Bessel functions, r, and r, are the inner and outer radii of the fin and B 


is given by: 





k is the thermal conductivity of the fin material and t 1s the half-thickness of the fin. The 


heat transfer coefficient from the fins is determined from: 


heDn _ (f/8)(Repn ee 1000)Pr 


(3.62) 
eo gos (hn) 


Nup, = 


where f is determined by: 
f = (0.79 In(Repn) — 1.64)” 


and the Reynolds number is: 





Repn = 


where the hydraulic diameter (D,,) is found from: 


D, = 228. 
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A is the flow area and P 1s the wetted perimeter as shown in Figure 3-8. The total heat 


transfer from the finned surface is: 
Q=h., Ad Tp — To): Ne (3.63) 


where A, is the surface area of the finned surface. This result is compared to the required 
heat flow per thermosyphon to determine if the finned condenser can be air-cooled. This 
completes the heat transfer circuit from the windings to the secondary coolant. For design 
purposes, the secondary coolant is assumed to have an entry temperature of 30 °C. The 


maximum allowable winding temperature is 150 °C. 
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Chapter 4: Experimental Setup 


4.1 Experiment Design 


The purpose of the experimental portion of this thesis is to determine if the 
theoretical results developed in Chapter 3 are accurate enough to be used for design 
purposes. With this objective in mind, the operating conditions and parameters are 
developed to specifically simulate the operating conditions that are expected to be 
encountered in ship propulsion motors. 

Due to size limitations of the laboratory facilities, a full scale experimental setup is 
impractical. Therefore, the experimental setup of a single RRT is designed to be a 
one-half scale model of the actual motor size. Similarity between the experiment and the 
actual motor is achieved by making the important non-dimensional parameters the same 
order of magnitude. Table 4-1 compares the important dimensions and non-dimensional 
parameters of the experimental setup and a typical destroyer size propulsion motor. 
Throughout this research the fill ratio is defined as the liquid fill volume divided by the 
evaporator volume. The length-to-diameter ratio is defined as the total thermosyphon 
length divided by the internal hydraulic diameter. 

All testing is completed by maintaining a constant rotational speed and increasing 
the power input to a predetermined value. The tests are run until a failure occurs or until 
a predetermined internal pressure of 50 psig is reached. Once failure occurs, the 


thermosyphon is cooled completely prior to re-starting it at a different speed setting. 
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Parameter Experiment 35,000 HP Motor 
Outer Radius (m) 
Inner Radius (m) 
Condenser Length (m) 
Evaporator Length (m) 
Rotational Speed (rad/sec) 
Film Reynolds No. = Re 


Archemedes No. = Ar, 


Ekman No. = Ek, 


Prandtl No. = Pr 


J vi/o 





Table 4-1: Experiment vs. Motor Parameters 


4.2 Measurements and Calculations 


In order to determine the heat transfer characteristics of the thermosyphon, the 


following quantities are measured: 


Cooling water inlet and outlet temperatures, T,, & T,,, [deg. C] 
Cooling water volume flow rate, Q,,,, [m’/sec. ] 
Speed of rotation, Q [rad./sec. ] 
Thermosyphon vapor temperature, T., [deg. C] 
Thermosyphon condenser wall temperature, T., [deg. C] 
Thermosyphon evaporator wall temperature, T., [deg. C]. 


Additional measurements taken are: 


Heater input voltage, V [ Volts] 
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- Heater input current, I [Amps] 
- Thermosyphon pressure, P,,, [kPa]. 


From these measurements the condenser and evaporator heat transfer coefficients and 
overall heat transferred by the thermosyphon are determined. The overall heat transfer 


rate comes from an enthalpy balance of the cooling water: 

Ome es elon — Vin) [Watts] (4.1) 
where C, is the specific heat of the cooling water. m is the mass flow rate of the cooling 
water which is: 


th = p20 « Qino [kg. / sec. ] (4.2). 


The cooling water volume flow rate is measured with a rotameter type flow meter. 
Density is determined from the saturated liquid thermodynamic Table at the average 


cooling water temperature. 


The heat transfer coefficients are calculated from the definition of a convective 


heat transfer coefficient: 
Q=h-A-(AT) [Watts] (4.3) 


where Q is determined from equation (4.1) above, A is the measured heat transfer area 
and AT is the difference between the thermosyphon wall and the saturation temperatures. 


This temperature difference is defined as: 


AT = (Tat — Tew) in the condenser 
AT = (Tew — Tsat) in the evaporator | 


Each of the above temperatures is measured in more than one way to obtain the 


final value. T_, is measured by two thermocouple probes. One probe is inserted into each 


sat 
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end of the test thermosyphon. These two measurements are averaged to obtain T,,. The 
pressure transducer output is used to compute T,,, independently and this value is 
compared to that obtained from the thermocouples as a sanity check. 

T.,, and T,, are each computed from the average of five wall temperature 
measurements. See Figure 4-1 for actual thermocouple locations on the test 
thermosyphons. All of the wall temperature thermocouples are located on the trailing 
edge of the thermosyphon for all tests. This location results in temperature measurements 
taken at the location of the greatest film thickness. Appendix A contains more detailed 


information on construction of the test thermosyphons. 
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Figure 4-1: Test Thermosyphon 
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T,,, and T,, are measured both in the fixed reference frame and on the rotor. Only 
the rotor measurements are used to calculate the heat transfer rate as they are more 
accurate than the stationary measurements. This is because the stationary measurement 
points are further from the test thermosyphon cooling jacket. The stationary 
measurements are uSed to set the cooling water flow rate since they are not subject to slip 
ring noise and can be read in real time. The flow rate is adjusted to maintain a cooling 


water temperature difference of approximately 10° C. 


4.3. The Test Apparatus 


A test stand was constructed to test thermosyphons in a rotating reference frame. 
The test stand is depicted in Figures 4-2 & 4-3. The main components of the test stand 
are the foundation, rotor, drive system, safety shield, cooling water system, heater and 
instrumentation. The important aspects of each will be described below. 

The foundation supports the rotor and drive system. It is designed to allow a rotor 
radius of 0.5 meter in keeping with the scale of the experiment. 

The rotor features rotating couplings on each end for the cooling water. There are 
also slip rings for power transfer to the heater and for connecting the rotating electronics 
assembly with the stationary one. A static balance of the rotor 1s sufficient as the 
rotational speed was limited to 500 rpm. 

The drive system consists of a DC motor which drives the rotor assembly through 
a V-belt with a 2.4:1 speed reduction. A variable transformer feeding through a bridge 
rectifier is used to vary the terminal voltage to the motor, achieving a continuously 


variable rotor speed. 
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Figure 4-2: Thermosyphon Test Stand 


The safety shield is shown in Figure 4-2. It 1s constructed of 1-1/2 inches of 
plywood with a 2"x4" frame surrounding the shield. The safety shield is sized to allow a 
rotor radius of 0.5 meter. Its strength is determined by assuming a 1.0 kg mass (weight of 
largest single rotating component) traveling at 5O m/s (maximum tip velocity) strikes the 
shield at mid-span between supports. Using the method described in [4.1] or [4.2] to 
calculate the depth of penetration provides a factor of safety of three or greater under 
these conservative assumptions. 

A schematic of the cooling water system is shown in Figure 4-4. Tap water is 
filtered and passed through a pressure regulator to reduce the pressure. A pre-heater 1s 


installed to allow positive control of the cooling water inlet temperature. A rotating 
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coupling passes the cooling water onto the rotor. Once there, it passes through a cooling 
jacket absorbing heat from the condenser section of the test thermosyphon. Another 

rotating coupling takes the water off the rotor where it passes through a second filter just 
prior to the flow meter. Upon exiting the flow meter, the cooling water is discarded into 


the city sewer drain. 
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Figure 4-3: Thermosyphon Test Stand 


The heater system consists of a regulated D.C. power supply connected to a 
resistance heater via carbon brushes and copper slip rings. Electric fans supply cooling air 
to the brushes and slip rings. The heater itself is a hollow cylinder of brass around which 
Nickel-Chromium resistance wire was wound. The resistance wire is potted inside two 


layers of thermally conductive ceramic. The heater is sized to produce 2.0 kW of power 
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at the maximum voltage rating of the power supply. It attaches to the outside of the test 
thermosyphon evaporator section via machine threads coated with heat-sink compound. 


The outside diameter of the heater 1s insulated to minimize heat loss to the surroundings. 
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Figure 4-4: Cooling Water System 


The instrumentation installed on the test stand consists of a flow meter (rotameter 
type), pressure transducer, optical speed sensor circuit, a pc-based 8-channel data 
acquisition system and a 16-channel thermocouple multiplexer/amplifier board. Figure 4-5 
shows a schematic of the instrumentation system. The multiplexer board is mounted on 


the rotor and is connected to the data acquisition system via copper slip rings and carbon 
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brushes. This arrangement sends three high-level (0-10 volt) analog signals across the slip 
rings: the multiplexed thermocouple signal, the pressure transducer signal and the cold 
junction compensation signal. The slip ring noise is canceled out by recording the voltage 


of a shorted input and subtracting this from the recorded input voltages. 
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Figure 4-5: Instrumentation Block Diagram 


Hand-held meters are used to read the input current and voltage as well as the 
shaft speed. These three measurements and the flow meter must be read manually, all 
temperature and pressure data are automatically recorded via the PC data acquisition 


system. 
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4.4 Accuracy and Error Analysis 


All measurements are imperfect attempts to determine an unknown quantity. 
Consequently, every measurement contains errors. These errors fall into two categories, 
systematic and random. Random errors are easily reduced by taking repeated readings. 
Systematic errors can only be eliminated by taking independent measurements using 
different techniques [4.3]. 

The systematic errors in the temperature measurements were reduced by 
calibrating the data acquisition system so that its output matched the reading of a 
calibrated glass bulb thermometer over the expected temperature range of 30-100° C. The 
flow meter was calibrated by measuring the volume of water passing through the meter in 
a fixed amount of time. 

The random errors were reduced by taking repeated readings and averaging. Data 
was collected over a 10 minute period for each data point. The manual readings were 
taken at one-minute intervals whereas the data acquisition system recorded all 
measurements at one-second intervals. 

Confidence limits on the temperature measurements are determined via standard 
Statistical methods. Hogg & Ledolter [4.4] shows that confidence limits are calculated 


from: 


nx ~ {s)-& a at ($)-£|=1-c (4.4) 


aoe 


JN 


(X) is 1 - a, where z(a/2) is the area under the standard normal distribution curve from 


which says the probability that the true mean (1) lies within z(a/2) of the sample mean 


minus infinity to o/2 and N 1s the number of samples taken. Equation (4.4) assumes that 
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the underlying distribution is normal with mean p and variance o*. For 95% confidence 
intervals a = 0.05 which makes z(a/2) = 1.96 [4.4]. 

A numerical example shows how this is applied to the measurements taken in the 
experiment. The selected temperature measurement is the worst-case encountered 1n the 
data. The measured inputs are the raw voltage, shorted voltage, cold junction voltage and 


slip ring voltage. For our example, these inputs are: 


Vraw = 5.40126 mV; s = 0.3826 mV 


V shorted = 1.60133 mV; s = 0.44998 mV 


Vije = 1.76722 V; s = 0.36407 V 


V stipring = 1.16734 V; s = 0.34671 V 


where N = 516 samples have been taken for each. Since N is large (greater than 30) the 
sample standard deviation (s) 1s assumed to be equal to the population standard deviation 


(o). Using equation (4.4) the 95% confidence intervals on these inputs are: 


aw = [5.368, 5.434] 
shorted = [1.563, 1.640] 
jc |= [1.736, 1.799] 
slipring — faliwhs.7 1197). 


ol ll 


The temperature of the cold junction is calculated from: 


Tg = 40.0. (V5. = ica) [deg. C] (4.5). 


In our example this 1s: 


eee 2 lo, zero | degC. 
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This temperature is entered into the National Bureau of Standards (NBS) thermocouple 
Tables to determine the cold junction voltage (y aare): The temperature at the 


thermocouple is calculated from: 


ie a Ves we Wanericd Bi Vcotdjana [mV] (4.6) 
This corrected voltage is entered into the NBS thermocouple equation to determine the 


measured temperature [4.5]: 
ie S06) ides. C 


It is emphasized that this + 3.3°C 1s a worst-case example. The vast majority of the data is 
within + 2.5°C or less. It is also worth noting that the majority of the uncertainty in the 
temperature measurement comes from the cold junction temperature. Since the heat 
transfer coefficients are calculated from temperature differences, rather than absolute 
temperatures, the uncertainty in the cold junction temperature has no effect on those 
calculations. 

A similar analysis carried out on the flow meter reading shows that the largest 
standard deviation (s=0.483) of the raw reading results in 95% confidence limits on the 
heat flow of +0.18 Watts/°C of cooling water temperature difference. For a typical 
temperature difference of 10 °C this gives confidence limits on the heat flow less than two 
watts. 

The above uncertainty can be carried through the calculation of the heat transfer 


coefficient or Nusselt number in non-dimensional form. The measured Nusselt number ts: 





= OL 
Nu = Si (4.7) 
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where k is the thermal conductivity, L and A are the characteristic length scale and heat 
transfer area, respectively. For the numerical example cited above, equation (4.7) 


evaluates to: 


Nu = 5737 + 1620 


which is a 95% confidence level of + 28% using the + 3.3°C temperature confidence 


limits. Again, this is the extreme case, the typical uncertainty is around + 20%. 
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41 US. Navy, Explosives Effects & Properties, (Confidential) NOLTR 65-218, 
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1987. 
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Chapter 5: Results and Conclusions 


5.1. Experimental Results 


Figures 5-1 through 5-15 at the end of this chapter summarize the experimental 
results. Figures 5-1 to 5-3 show in dimensional form, the effects of effective gravity, 
working fluid and length to diameter ratio in the condenser. The overall trend in these 
graphs shows that heat flux is proportional to temperature difference. 


Recall that the definition of a convective heat transfer coefficient 1s: 


mace : 
baa (5.1). 


Figure 5-1 shows that the heat transfer coefficient increases with effective gravity as was 
predicted in Chapter 3. This is observed by noticing that for any given temperature 
difference, the measured heat flux increases with relative gravity. 

The effect of working fluid 1s depicted in Figure 5-2. For a given temperature 
difference, the heat flux is higher with methanol than with ethanol. This is a consequence 
of the higher heat of vaporization of methanol. This is also predicted from the results of 
Chapter 3. 

Figure 5-3 shows the effect of the L/D ratio on the condenser. The only 
discernible difference between the three L/D's tested 1s that the smaller one (L/D = 30) 
exhibits a large amount of scatter in the data whereas data from the larger L/D ratio 
thermosyphons are more consistent. This is probably attributable to two-dimensional 


effects which cannot be predicted with the one-dimensional model used in Chapter 3. 
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Figures 5-4 to 5-7 show the same data in non-dimensional form. These results 
confirm the relationship predicted by equation (3.26) for the condenser heat transfer 
coefficient. Figures 5-4 and 5-5 show that the effects of both the effective gravity and the 
working fluid are accounted for by equation (3.26). The experimental results appear to be 
about 20% higher than the theoretical predictions. However, Figures 5-6 and 5-7 show 
that this 1s largely due to the use of the smaller (L/D = 30) thermosyphon where the 
two-dimensional effects come into play. 

Figure 5-8 shows the effect of the fill ratio on the maximum heat flux. 
Thermosyphons were tested at fill ratios from 11% to 85%. Recall that the fill ratio 1s 
defined as the ratio of liquid volume to evaporator volume. At both extremes, failure 
occurred at very low heat fluxes and showed little improvement with increased effective 
gravity. In the intermediate range, there was significant improvement in the maximum 
heat flux with increased effective gravity. There appears to be an optimum fill ratio in the 
vicinity of 35-45%. At the lower fill ratios, the early failure is caused by evaporator 
dry-out. At the higher fill ratios, the early failure is caused by condenser flooding. This 1s 
consistent with previous results on gravitational thermosyphons published in references 
[5.1] and [5.2]. The filling method employed did not allow precise control over the fill 
ratio so this portion of the results is considered very approximate. See Appendix A for 
more discussion of the filling method. 

The evaporator data is plotted in Figures 5-9 to 5-14. The dimensional data is 
presented in Figures 5-9 to 5-11. Because of the large amount of scatter in this data, the 


regression line on these figures is of little value and not much can be said about Figures 
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5-9 and 5-10. Figure 5-11 does exhibit the same behavior as that of Figure 5-3 for the 
condenser. That is, the scatter in the data is reduced at the higher L/D ratios. Eliminating 
the L/D = 30 data from Figures 5-9 and 5-10 will result in a relationship between heat flux 
and temperature difference similar to that found in Figures 5-1 and 5-2 for the condenser. 
This is shown in Figure 5-16. 

The nondimensionalized data in Figures 5-12 to 5-14 clearly shows the evaporator 
heat transfer is a combination of pool boiling as predicted by equation (3.44) and film 
evaporation which is predicted by equation (3.33). The data shows that equation (3.33) is 
fairly accurate for Nusselt numbers above about 2500. When equation (3.33) predicts a 
Nusselt number below 2500, the actual Nusselt number remains fairly constant at that 
value. This suggests the following relationship to predict the evaporator heat transfer 


coefficient: 


0.82 
a Pr | Rel-22 _ 3a Netley" one? 
4Ja L 


Pr 4/3 
Nut, = Set = max ( {ve ) (5.2). 
l 


2500 


Therefore, equations (4.8) and (3.26) should be used to predict the evaporator and 
condenser heat transfer coefficients, respectively. 

In order to examine the extent of the Coriolis effect, one set of tests was run with 
the direction of rotation reversed so that the wall temperature thermocouples were located 
on the leading vice the trailing edge. As described in Chapter 3, the Coriolis effect causes 
the film to be thinner on the leading edge. This is demonstrated in Figure 5-15 by a 


smaller temperature difference for a given level of heat flux. Figure 5-15 also shows that 


99 


this effect 1s rather small as expected. It is also interesting to note that dry-out failure was 
observed at a lower heat flux when the thermocouples were located on the leading edge. 
This phenomenon is also predicted in Chapter 3. All data gathered in the experimental 
portion of this thesis are located on the enclosed floppy disks in Lotus 123 Release 4 


format. A summary of the data is tabulated in Appendix B. 


5.2 Preliminary Design Tool 

A computerized preliminary design tool has been developed for sizing multi-disk 
axial gap propulsion motors to a specific ship application. This allows the ship designer to 
easily conduct trade-off studies to determine the best propulsion plant configuration for a 
particular set of ship requirements. Figure 5-17 shows the input and output sections of the 
design tool. The heat transfer section is summarized in Figure 5-18. Stator cooling is via 
indirect water and rotor cooling is via radial rotating thermosyphons. Air cooling of the 
thermosyphon condensers is evaluated both with and without fins for comparison with 
water cooling of the condensers. Air gap heat transfer is also considered in the heat 
transfer model. Other heat transfer paths are negligible by comparison and are ignored. 
The equations developed in Chapters 2 and 3 constitute the core of the design tool. 

The design tool was wnitten as a spreadsheet in the Lotus 123 release 4 program. 
A spreadsheet format was chosen because it allows the user to change any input and 
immediately see the effect on any of the outputs. The spreadsheet has four pages. Figures 
5-17 and 5-18 constitute pages one and two, respectively. Page three contains various 


property data and page four is a record of changes log. The shaded blocks in Figure 5-17 
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constitute the design inputs. Additional inputs are the empirical boiling coefficients 
located in the shaded area in Figure 5-18 and the various material properties located in 
shaded regions of page three of the spreadsheet (not pictured). 

Additionally, the user must use the "solver" feature of the spreadsheet program to 
solve the two implicit equations in the spreadsheet. The first comes from equation (2.38) 
in the structural model. The key sequence AL7-R-A-S brings up the solver dialog box. 
Cell A:B49 1s selected as the adjustable cell and cell A:B50 is the constraint cell. There 1s 
no optimal cell. Select enter and the spreadsheet will solve the implicit equation using a 
Newton-Raphson technique. The second implicit loop comes about from the heat transfer 
model. The thermosyphon design section uses the total heat transfer required as an input 
to calculate the winding temperature whereas the air gap heat transfer section uses the 
winding temperature and the incoming air temperature to calculate the total amount of 
heat transfer obtained. Since other heat transfer paths are neglected, conservation of 
energy requires the total heat transferred via these two paths to equal the total losses 
generated. This sequence of equations creates an implicit loop in the spreadsheet. This 
loop is also solved with the "solver" function by specifying cell B:E8 as the adjustable cell 
and B:E12 as the constraint cell. With all of the inputs specified and the two implicit 
equations solved, the design is complete at this level of detail. A copy of the design tool 1s 


included on the disks attached to this thesis. 


101 





Synchronous Machine Design Program 


Date: 5-8-95 Axial Fleld Motor 
Inputs are shaded 
Design Parameters Variables of interest 
erminal Voltage 7530 “Volts eaf 2.00 per unit 
Rated Speed 164 mm xd 1.35 per unit 
Rated Frequency 188.5 rad/sec Jas 7.4E+06 A/m‘42 
Rated P.F. angle 0.451 = rad Jfs 7.7E+06 A/m‘*2 
Outer radius, Ro 1.08 m volume 5.14 m3 
Inner radius, Ri 0.65 m weight 20.3 LT 
Air gap, g 0.01 ™m Efficiency 96.2% 
No. of air gaps, ng 10 Traction 24.61 psi 
Pitch, p 0.889 Power 2.6E +07 WwW 
No. of slots/pole - phase, q 3 Power 35355 HP 
tator space factor, lars 0.7 Stack Length 0.846 m 
Stator slot depth, hs 0.027 m 
tator slot factor, SSF 0.8 
Lams at outer radius = 0.42 m 
Core lamination thickness 0.0005 m Stator 
No. rotor slots/pole, qr 4 # of pole pairs, p 11 
Rotor slot depth, hr 0.022 m No. of series stators, Ss 10 
Rotor space factor, lamr 0.7 Pitch factor, kp 0.985 
Rotor slot factor, RSF 0.8 Breadth factor, kb 0.960 
Lamr at outer radius = 0.42 m Voltage, Va 102 Volts/tum 
Stator electric loading, Ks 2.0E+05 A/m Current, la 95295 Amp-tums 
Peak flux density in iron, Bs 18 TT Rated Apparant Power 2.93E+07 V-A 
Peak air gap flux density, Bgap ae 6|OUdT Base Impedance, Zb 0.00108 Ohms/tum*2 
Max. cooling air velocity, Vs 25 ms Slot current density, Jas 7.4E+06 A/m‘2 
Cooling Air entry temp., To 30 deg.C Cu current density, Ja 9.3E+06 A/lm2 
Tsat 100 deg. C Tum Length, It 1.30 meters 
Thermosyphon Diameter, D 0.007 m Resistance, Ra 2.14E-05 Ohms/tum’2 0.0199 per unit 
Fill Ratio 0.2 Flux per pole, phi 0.08133 Webers 
Insulation Thickness 0.005 m Back iron depth,tb 0.105 meters 
Insulation thermal Cond. 2 Wim-K Slot width, ws 0.014 meters 
Fin Spacing 0.01 m 
Fin Thickness 0.001 m Rotor 
Thermal Conductivity of Alum. 111 W/mK No. of senes rotors, Sr 10 
Tip speed, u 18.5 m/s 
Structural Design Rotor breadth factor 0.958 
Shaft Outer Diameter 0.509 m Radus ratio, a 0.60 
Shaft Inner Diameter 0.408 m excitation, eaf 2.00 per unit 
Shaft Cross Section 0.07339 m**2 tum length, Itf 1.36 meters 
Torque per Disk 307026 N-m field resistance 9.1E-06 Ohms/tum’2 0.0085 per unit 
Structural Outer Radius, r_oo 1.39 ™m N/L field 121563 Amp-tums 
Structural Inner Radius, r_ii 0.46 m 
Rotor Disk Thickness 0.0023 m Full load field 243402 Amp-tums 
Stator Disk Thickness 0.0018 m 
Drum Inner Diameter 0.904 mM Cu current density, Jf 9.7E+06 A/m*2 
Case Outer Diameter 2.78 =m Slot current density, Jfs 7.7E+06 A/m*2 
Solver Constraint Cell 0 Electnc loading, Kf 1.7E+05 A/m 
Welght and Volume Slot width, wr 0.033 meters 
Stack length, L 0.846 m Reactances 
Air gap volume 0.234 m‘3 Leakage 5.77E-05 Ohms/tum’2 0.05 per unit 
Back iron vol. 0.491 m3 Xaao 9.31E-04 Ohms/tum*’2 0.87 per unit 
Stator iron vol. 0.395 m‘3 D-axis, xd 1.45E-03 Ohms/tum’2 1.35 per unit 
Stator Winding volume 0562. m*3 Losses 
Rotor Winding volume 0.426 m3 Stator copper 583325 Watts 0.0199 per unit 
Rotor steel vol. 0.677 m3 Windage & frict. 35 Watts 0.0000 per unit 
Case vol. 0.019 m‘3 Core loss 24741 Watts 0.0008 per unit 
Shaft volume 0.062 m3 Field copper 540899 Watts 0.0185 per unit 
Machine volume 5.145 m3 Total Losses 1149000 Watts 0.0392 per unit 
Rotor weight 9077 kg Efficiency 96.2% 
Stator weight 7543 kg Heat Transfer Summary 
Back iron weight 3831 kg Stator air-gap heat loss 136708 Watts 
Case weight 146 kg Rotor air-gap heat loss 302870 Watts 
Machine weight 20597 kg Rotor t'syphon heat loss 238029 Watts 
Machine weight 20.3. LT Stator water heat loss 471359 Watts 


Figure 5-17: Design Tool Input and Output Sections 
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Figure 5-18: Design Tool Heat Transfer Section 
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5.2.1 Design Example 


The inputs of Figure 5-17 are chosen to produce a small diameter 35,000 HP 
motor suitable for a podded propulsion destroyer application. It is a 10 air gap or 5-disk 
design. Figure 5-19 1s a scaled preliminary drawing of this design. Figures 2-5 and 2-6 
also depict this specific design. Table 5-1 compares several other existing or proposed 
motor technologies at the same power and speed rating. Note the significant improvement 
in both volume and weight over the other technologies. The numerical values for the 
other technologies were obtained from [5.3]. 

In addition to the values obtained from [5.3], the permanent magnet single disk 
axial gap motor was simulated with the design tool. It is interesting to note that this 
motor 1s actually only a 30,000 HP motor vice the 35,000 HP advertised. Comparison of 
the starred (*) values calculated with the design tool and non-starred values for this motor 
show that the design tool is giving reasonable results. The weight difference is due to the 
difference between permanent magnet and wound field synchronous technologies. The 
length difference is accounted for by the fact that the design tool does not include a length 
allowance for journal or thrust bearings which are included in the other data. This length 
difference also accounts for the volume difference. 

Comparison of the starred values in the last two columns of Table 5-1 show the 
advantages of the multi-disk technology over the next most recent proposed motor design. 
The multi-disk axial design provides a 29% weight reduction, a 40% volume reduction 
and nearly a 50% reduction in radius over the single-disk design while providing 16% 


more power. 
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Conventional Conventional PM axial gap Multi-disk 
air-cooled = water-cooled motor Axial motor 


Cooling method air/air Water/water water/air —_ water/t'syphon 
(Arm / Field) 


Armature current 6.2 Ue 7.4 
density (A/mm’) 


Field current 1#s ey 
density (A/mm’) 


Airgap flux density 1.2 0.84 2 
(T) 

Electromagnetic lS8eZ 

shear stress (KPa) 


) 98.5 
eal 
L204" 
22 / 28.1* 
22.7 / 8.4* 
Table 5-1: Comparison of Motor Technologies 





5.3 Conclusions and Recommendations for Future Research 


Table 5-1 shows the significant volume and weight advantages of the multi-disk 
design developed in the course of this research. Of particular interest to the ship designer 
is the nearly 50% reduction in the motor diameter with the multi-disk geometry. This 
technology clearly shows promise and more research should be conducted in this area. 

The computer model developed herein is only the first step in evaluating this new motor 
geometry. It will however, allow ship trade-off studies to be completed to determine if the 
ship impacts predicted in the introduction are favorable enough to make electric drive cost 
competitive. Such trade-off studies are the first recommendation for future research. 

The motor design tool developed in the course of this research is very well suited 


to being incorporated into an optimization program. This feature allows future 
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Figure 5-19: Motor Section Views 
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researchers to explore the entire design space of these unique motors, eventually 
developing an optimum geometry for ship propulsion motors. Once the preliminary design 
is optimized, a detailed design sequence should be undertaken. The detailed design should 
include a complete structural analysis and a more detailed electromagnetic analysis which 
includes harmonic effects, fringing fields and end turn effects. Once these analyses are 
completed, a prototype multi-disk motor should be built and tested. 

This thesis has presented theoretical predictions with experimental verification of 
the heat transfer coefficients found in radial rotating thermosyphons. While this research 
is considered successful because of the agreement between theory and experiment, the 
application of thermosyphons to rotor cooling cannot yet be recommended without 
reservations. Specifically, the experimental portion of this research revealed that there are 
serious stability problems with the RRT. 

It was noted during testing that the heat flux at failure was not path independent. 
Failure could be induced at different power levels for the same speed setting depending on 
how the operating point was approached (i.e. holding speed constant and increasing 
power or changing speed at a constant power level). Because of this hysteretic behavior, 
a prediction for the maximum heat flux could not be developed. Such a prediction 1s a 
necessity if optimal designs are to be generated. Once failure occurred, the thermosyphon 
had to be cooled to room temperature prior to restarting. 

Transient testing was also not part of this research. Transient operation is an 
important part of the operating profile of a ship propulsion motor. It is recommended that 


any future research into radial rotating thermosyphons include transient effects. 
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Other areas which should be considered in future research are a more detailed 
study of the effects of fill ratio on the maximum heat flux, larger length to diameter ratios 
and other working fluids. It is believed that the theoretical predictions presented herein 
will exhibit good agreement with experiment for L/D ratios of 50 or greater. The motor 
design tool currently uses water as the working fluid. If the problems of freezing can be 
overcome, water will be a superior working fluid to the alcohols because of its very high 
heat of vaporization. 

The shape of the thermosyphon cross-section should also be evaluated. All testing 
in the current research was carried out on circular cross-sections. However, as shown 1n 
Figure 3-6 a rectangular cross-section would be more effective for conducting the heat 
into the evaporator. Depending on the specific section chosen, the Coriolis effect may 
become more significant under operating conditions encountered in ship propulsion 


motors. 
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Condenser Temperature Difference vs. Heat Flux 
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Figure 5-2 
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Condenser Temperature Difference vs. Heat Flux 
Effect of L/D Ratio 


= 
© 
Nn 


~-_~ 
N 
= 
— 
= 
x 
5 
Le 
al 
O 
® 
a0 


QO A OOD 


Cm O 


Temperature Difference (deg. C) 


















= 
XN prediction 
J . _99 
-p 4 4 10g 
- = 
ri v _17g 
Oo 2 ° _23g 
E i _28g 
D 2 : _32g 
5 : 46 
S _*99 
_99g 






iets gineel: yf (or 






10° | 
10° 2 3 4 5 


Nu, = 3/4*(Pri(4*Ja))° 8*4((X+b,)43-x 43 \(sqrt( v/a) /3}° 82 
Figure 5-4 
111 









Condenser Predicted vs. Measured Nusselt No. 
Effect of fluid 
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Condenser Predicted vs. Measured Nusselt No. 
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Condenser Temperature Difference vs. Heat Flux 
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Appendix A: Thermosyphon Construction 


This appendix describes the techniques and lessons learned in construction and 
operation of the test thermosyphons. The primary references used were Dunn & Reay 
[A-1] and Chisolm [A-2]. Both of these references contain fairly detailed descriptions of 
construction techniques, material selection and compatibility information and a list of 
commercial manufacturers of heat pipes. Dunn & Reay also contains a chapter on the 
practical aspects of heat pipe design which was very useful in designing the experiment. It 
is the intent of this appendix to augment rather than repeat information found in those 
references. 

Selection of working fluids was based primarily on the expected operating 
temperatures of the thermosyphon and the fluid properties. Water, Methanol and Ethanol 
were chosen as test fluids. Water was later abandoned because of the difficulty in its use. 
There were two problems with water: 1) removal of dissolved gasses was extremely 
difficult and; 2) the low vapor pressure of water when at room temperature placed very 
rigorous requirements on the various seals of the thermosyphon to prevent air from 
leaking into It. 

The container material was initially selected based upon compatibility with the 
working fluid. Dunn & Reay recommended copper for use with these working fluids. 
Because of its very high thermal conductivity, copper caused problems obtaining accurate 
temperature measurements. It is believed that this was caused by axial conduction of heat 


in the container walls. The thermosyphons from which the data was gathered were 
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constructed from yellow brass and type 304 stainless steel. Dunn & Reay reported that 
some researchers had experienced problems with stainless steel and the alcohols, however 
this author experienced no problems with stainless steel. Brass was a good compromise 
material since its thermal conductivity is one fourth that of copper but it is much easier to 
solder than stainless steel. 

For this test program, the ability to empty and re-fill the thermosyphons was 
desired. Consequently, the filling valve was left on the thermosyphons during testing to 
seal out the atmosphere. This worked quite well. Brass plug valves with Viton O-rings 
were used quite successfully. New valves were disassembled and cleaned with acetone. 
The seals were then coated with silicone vacuum grease prior to re-assembly. 

Both of the references cited above recommend welding or soldering for sealing 
heat pipes. All mechanical joints were sealed by soldering with Stay-Brite silver bearing 
solder. This solder's melting point is lower than normal silver solders, but is higher than 
lead-tin soft solder. It also has much greater strength characteristics than soft solder at the 
operating temperatures encountered in this research. 

The fill valve, pressure transducer and the internal thermocouple probes were 
attached to the thermosyphon via mechanical connections to allow for removal and 
cleaning between uses. Swagelock brand connectors initially were used for these 
mechanical connections. Unfortunately, after assembling and disassembling these 
connectors 2-3 times they became susceptible to leaks. The thermosyphon design was 
then altered to accommodate Cajon VCO vacuum fittings. These fittings worked 


extremely well. In addition to being leak free and more sturdy, they allowed larger access 


119 





holes into the thermosyphon for brush cleaning the interior. Figure 4-1 depicts this 
version of the test thermosyphon. 

The cleaning and filling procedure was adapted from several described [A-1] and 
[A-2]. The need for cleanliness when assembling thermosyphons cannot be 
overemphasized. Removal of all oils, greases and oxides from the container is essential to 
maintaining a vacuum inside the thermosyphon after disconnecting it from the vacuum 


system. The cleaning procedure used in this research was as follows: 


1. Soak all components in a base bath overnight. Base bath formula is one cup of 
Potassium-hydroxide chips, one cup of distilled water and two gallons of 
methanol. 

2. Passivate components by soaking for 15 minutes in a 15% nitric, 5% hydrochloric 

acid bath. 

Rinse in distilled water. 

After assembly, rinse inside of thermosyphon with trichloroethane. 

Soak for 2 hours in trichloroethane. 

Rinse out with acetone. Clean out with brush soaked in acetone. 


Soak for 2 hours in acetone. 


OCS = = 


Rinse with distilled water. 


Steps 4-8 were repeated each time a thermosyphon was re-filled. 
The filling procedure 1s also critical in maintaining a vacuum long enough to 


complete testing. The filling procedure was: 


1. Maintain vacuum in filling ng at all times. 
Evacuate a new thermosyphon for a minimum of two days. 
Close thermosyphon fill valve and vacuum pump isolation valve. 


Remove and weigh evacuated thermosyphon. 


waco a 


Re-install thermosyphon, open vacuum pump isolation valve. 
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6. When ultimate vacuum is achieved, open thermosyphon fill valve. 

7. Evacuate thermosyphon until there is no difference in ultimate vacuum with the fill 
valve open or closed. 

8. Isolate vacuum pump from system by closing valve. 

9. Open isolation valve to liquid burette. 

10. Meter desired amount of fluid into fill ng. 

11. Heat walls of fill mg to evaporate liquid off walls. 

12. Close burette isolation valve and thermosyphon fill valve. 


13. Remove and re-weigh thermosyphon to determine exact amount of liquid fill. 


The design of the fill rig was essentially identical to those described in [A-2]. 
Figure A-1 depicts the configuration of the fill rig. The problem encountered with the fill 
rig was the inability to install a precisely measured amount of liquid into the test 
thermosyphon. Some of the fluid coats the walls of the fill ng, thus the volume 
measurement from the burette is only approximate. The most accurate measurement of 
the liquid in the thermosyphon comes from the weight difference. However, this is done 
after the fact so the only way to remedy an incorrect fill is to empty the thermosyphon and 
start over. 

Type E thermocouples were used throughout the experiment. Initially, silicon 
diaphragm pressure transducers were used as the manufacturer indicated they would be 
acceptable for use with distilled water and the alcohols. However, the distilled water 
caused the silicon to rupture in a very short time, creating vacuum leaks in the 
thermosyphons. Pressure transducers with all stainless steel wetted parts were substituted 


with satisfactory results. The specific model used (Omega Engineering PX-120), did have 
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Figure A-1: Thermosyphon Fill Rig 


a tendency to take a permanent set altering its zero reading significantly after exposure to 
a vacuum or vapor pressure for a few days. This created difficulty in correlating the 
pressure and temperature measurements until the thermosyphon was emptied and the zero 
on the transducer could be re-checked. Once the proper zero was entered into the data 
sheets, good correlation between the measured saturation pressure and measured 


saturation temperature was obtained. 
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A-] Dunn, P.D. & Reay, D.A., Heat Pipes, 4th ed., Pergamon, 1994. 


A-2 Chisolm, D. et al., "Manufacture of Heat Pipes," Analysis Design & Manufacture 
of Heat Pipes, National Engineering Laboratory Course Papers, 19-21 June, 1973. 
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Appendix B: Data Summary 


gravity 
(g's) 
1] 

11 

1] 

1] 

1] 

11 

1] 

1] 


q_cond 
(W/m“2) 
32,007 
Dy o20 
(WP 
91,758 
148,160 
200,038 
256,708 
329,940 
55,081 
92,791 
150,845 
198,899 
BSS LOTS 
299,693 
52,666 
86,757 
144,089 
181,808 
249,425 
284,350 
47,963 
85,860 
124,299 
176,307 
239,314 
276,588 
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Methanol, L/D=53, F/R=0.40 


t_sat 
(deg. C) 
45 
60 
68 
80 
94 
110 
117 
148 
58 
73 
89 
101 
118 
133 
52 
68 
81 
95 
110 
124 
53 
67 
a7 
92 
105 
HK? 


aelatac 
(deg. C) 
4 
8 
10 
13 
21 
30 
42 
60 
6 
i 
20 
28 
37 
48 
6 
10 
he 
23 
30 
39 








Methanol, L/D=53, F/R=0.44 


pipe no gravity q_cond t_sat dele Gem € 
(g's) (W/m’2) (deg. C) (deg. C) (deg. C) 

1] 11 22,862 42 3 6 
11 JJ 36,674 49 6 7 
1] 1] 59,572 56 10 a 
te 11 97,506 68 14 10 
1] 11 114,505 77 18 eZ 
11 11 138,489 86 23 13 
1] 1] 182,010 oy 26 18 
te 16 39,148 49 6 6 
11 16 87,410 64 12 10 
a 16 109,098 75 17 2 
ih 16 166,008 86 24 13 
I 16 190,171 96 38 16 
1] 28 18,978 46 5 14 
1] 28 67,654 64 11] 11 
11 28 95,935 71 15 13 
11 28 2770 78 19 14 
11 28 172,667 86 24 16 
11 28 197,517 104 32 18 
11 28 258,194 113 38 21 
11 47 Z5e2N9 42 * 13 
1] 47 62,168 60 9 10 
1] 47 86,669 68 12 14 
11 47 135,687 81 17 17 
Ie 47 352351 94 24 19 
1] 47 247,157 106 3] zee 
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Methanol, L/D=53, F/R=0.36 





pipe_no gravity q_cond t_sat del tc del te 
(g's) (W/m”2) (deg. C) (deg. C) (deg. C) 
12 9 22,815 46 4 4 
12 9 38,703 54 6 6 
12 9 60,961 59 10 8 
12 o 86,102 7) 15 9 
12 9 115,705 8] 20 10 

12 16 23,821 45 4 

12 16 38,628 8 

12 16 59,865 58 9 8 
12 16 87,445 66 13 9 
FZ 16 112,197 74 17 10 
12 16 134,965 85 Zi 11 
12 16 184,152 96 26 12 
12 23 39,035 48 5 6 
12 23 62,515 58 9 g 
12 23 93,404 70 14 10 
12 23 138,398 82 20 10 
12 23 192,092 99 Zi 1] 
12 BB 235,138 107 31 13 
(2 23 256,450 LZ7 39 13 
12 46 40,787 0 i 18 
| 12 46 75,159 61 10 15 
12 46 124,191 15 Ld 14 
1 46 170,216 90 23 1S 
He 46 245 103 29 18 
12 46 282,657 12 36 20 
12 46 309,23 1 123 oo 20 
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Methanol, L/D=30, F/R=0.33 


pipe no gravity q_ cond t_sat aa te acinese 
(g's) (W/m*2) (deg. C) (deg. C) (deg. C) 

13 18 24,249 46 2 7 
13 18 34,309 54 5 6 
13 18 57,347 59 6 9 
13 18 80,193 70 9 1] 
13 18 98,207 oo 1] 12 
13 18 121,136 87 15 2 
13 18 155,991 93 17 1] 
13 18 212,319 107 24 9 
13 32 49,304 60 5 13 
13 32 68,606 63 6 LS 
13 32 89,075 70 9 16 
13 32 148,820 86 14 12 
13 32 188,878 9] 18 11 
13 45 53,693 56 4 19 
13 45 71,700 60 6 20 
13 45 101,958 71 9 Z| 
13 45 148,238 80 12 35 
13 45 148,512 79 12 18 
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Methanol, L/D=30, F/R=0.25 


pipe _no gravity q_cond t_sat Gelst ic del t_e 
(g's) (W/m‘’2) (deg. C) (deg. C) (deg. C) 

14 1] 38,278 58, i) 14 
14 1] 54,692 64 8 IVF 
14 1] 70,951 67 1] 18 
14 1] 106,335 75 14 16 
14 11 120,667 83 17 18 
14 1 1527750 94 23 20 
14 16 52,328 60 8 17 
14 16 71,524 66 10 7 
14 16 106,575 76 14 16 
14 16 138,256 88 18 24 
14 28 72,238 62 ql 21 
14 28 95,470 69 10 oa 
14 28 134,265 80 15 26 
14 28 168,817 97 32 29 
14 47 61,073 66 5 23 
14 47 87,731 68 9 26 
14 47 141,369 76 he 34 


Ethanol, L/D=53, F/R=0.85 


pipe no gravity q_cond t_sat del t_c del te 
(g's) (W/m‘’2) (deg. C) (deg. C) (deg. C) 

15 8 137 44 0 28 

les 30 Deon 50 12 19 

Ms) 53 9,172 47 1] 18 

NS: Wa 1,470 38 5 10 

15 75 9,231 a6 16 28 

15 75 23,798 54 i 30 

15 dS 49,059 58 8 32 
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Ethanol, L/D=30, F/R=0.20 


pipe_no gravity q_cond t_sat del tc del t_e 
(g's) (W/m’2) (deg. C) (deg. C) (deg. C) 

ey 6 13,020 42 6 18 
ly 6 259 51 6 Zl 
17 6 39,440 58 9 Zl 
17 6 50,446 58 10 oz 
17 6 61,131 68 14 22 
17 10 15,055 43 6 14 
17 10 27,767 49 6 25 
17 10 34,192 55 8 23 
17 10 42,715 56 9 Ze 
17 10 54,220 65 13 23 
17 15 9,243 45 8 15 
17 I 26,025 48 6 2S 
hy 15 23,973 Sy). 6 25 
17 15 43,086 56 8 24 
Ne? 15 59,392 64 11 23 
17 59 12,383 43 5 26 
17 59 21,796 46 4 29 
17 59 40,553 56 5 30 
Ie 59 54,573 60 8 31 
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Ethanol, L/D=30, F/R=0.28 


pipe _no gravity q_cond t_sat Getetc Gclarae 
(g's) (W/m‘’2) (deg. C) (deg. C) (deg. C) 

18 6 22 > 6 10 
18 6 356070 60 9 2 
18 6 60,917 70 14 13 
18 6 54,351 78 18 ie 
18 6 99,908 90 24 2 
18 6 120,787 99 31 I 
18 6 189,956 120 47 1] 
18 6 232,413 143 66 1] 
18 10 19,860 52 7 20 
18 10 41,784 59 10 13 
18 10 66,726 72 eS 13 
18 10 120,744 9] Ze 14 
18 10 186,091 118 43 eS 
18 10 251,560 142 61 10 
18 15 19,929 a3 8 18 
18 IS 58,328 64 10 14 
18 15 90,256 80 18 Is 
18 15 145,827 99 30 14 
18 15 ZO SOS ey 42 13 
18 is 252,768 134 54 10 
18 59 35,790 59 6 24 
18 59 T3715 70 10 24 
18 59 125,028 84 18 4] 
18 59 184,305 96 25 19 
18 59 230,313 113 34 16 
18 59 290,267 124 42 16 
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Methanol, L/D=62, F/R=0.29 


pipe _no gravity q_cond t_sat Geluine delle 
(g's) (W/m‘2) (deg. C) (deg. C) (deg. C) 

19 16 77,098 63 13 20 
19 16 93,262 75 17 21 
19 28 55,160 54 8 12 
19 28 79,707 65 12 16 
19 28 99,893 1 16 19 
19 28 131,986 81 20 oo 
19 28 154,973 9] 29 24 
19 47 46,676 58 8 1] 
19 47 77,857 61 10 15 
19 47 100,594 70 14 16 
19 47 131,421 74 16 2| 
19 47 168,541 84 22 24 
19 47 214,438 92 27 26 
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pipe no 


12a 
12a 
12a 
12a 
12a 
12a 
12a 
12a 
12a 
12a 
12a 
12a 
12a 
12a 
12a 
12a 
12a 
12a 
12a 
12a 
12a 
12a 
12a 
12a 


Methanol, L/D=53, F/R=0.36 Reverse Rotation 


gravity 


(g's) 


q_cond 
(W/m“2) 
27,545 
32,878 
61,189 
21,500 
40,658 
60,663 
84,704 
18,764 
33,46] 
54,775 
82,520 
105,826 
35,006 
64,430 
92,241 
122,350 
183,890 
47,821 
74,047 
ieee 
182,529 
229,534 
285,908 
288.392 
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t_sat 
(deg. C) 
48 
6] 
70 
49 
51 
6] 
74 
43 
47 
56 
68 
79 
46 
60 
70 
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99 
50 
59 
ie 
90 
105 
119 
124 


gdelitec 
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